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Prologue 

This lecture courso0 is intended to fill the gap between graduate courses on 
quantum field theory and specialized reviews or forefront-research articles 
on functional renormalization group approaches to quantum field theory and 
gauge theories. 

These lecture notes are meant for advanced students who want to get ac- 
quainted with modern renormalization group (RG) methods as well as func- 
tional approaches to quantum gauge theories. In the first lecture, the func- 
tional renormalization group is introduced with a focus on the flow equation 
for the effective average action. The second lecture is devoted to a discussion 
of flow equations and symmetries in general, and flow equations and gauge 
symmetries in particular. The third lecture deals with the flow equation in the 
background formalism which is particularly convenient for analytical compu- 
tations of truncated flows. The fourth lecture concentrates on the transition 
from microscopic to macroscopic degrees of freedom; even though this is dis- 
cussed here in the language and the context of QCD, the developed formalism 
is much more general and will be useful also for other systems. Sections which 
have an asterisk * in the section title contain more advanced material and 
may be skipped during a first reading. 

This is not a review. I apologize for many omissions of further inter- 
esting and important aspects of this field (and their corresponding refer- 
ences). General reviews and more complete reference lists can be found in 
[ll[2l[3l[4l[5l[6l[7l[8]. A guide to more specialized literature is given in the 
Further- Reading subsections at the end of some sections. 



1 Introduction 

Quantum and statistical field theory are two sides of the same medal, repre- 
senting the fundament on which modern physics is built. Both branches have 
been molded by the concept of the renormalization group. The renormalization 

^ Lectures held at the 2006 ECT* School "Renormalization Group and Effective 
Field Theory Approaches to Many-Body Systems" , Trento, Italy. 
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group deals with the physics of scales. A central theme is the understanding 
of the macroscopic physics at long distances or low momenta in terms of the 
fundamental microscopic interactions. Bridging this gap from micro to macro 
scales requires a thorough understanding of quantum or statistical fluctuations 
on all the scales in between. 

Field theories with gauge symmetry are of central importance, in partic- 
ular, since all elementary particle-physics interactions are described by gauge 
theories. Moreover, nonabelian gauge theories - and most prominently quan- 
tum chromodynamics (QCD) - are in many respects paradigmatic, since they 
exhibit numerous features that are encountered in various field theoretical 
systems both in particle physics as well as condensed-matter systems. During 
the transition from micro to macro scales, these theories turn from weak to 
strong coupling, the relevant degrees of freedom are changed, the realization 
of the fundamental symmetries is different on the various scales, and the phase 
diagram is expected to have a rich structure, being formed by different and 
competing collective phenomena. 

A profound understanding of gauge theories thus requires not just one but 
a whole toolbox of field theoretical methods. In addition to analytical per- 
turbative methods for weak coupling and numerical lattice gauge theory for 
arbitrary couplings, functional methods begin to bridge the gap, since they are 
not restricted to weak coupling and can still largely be treated analytically. 
Functional methods aim at the computation of generating functionals of cor- 
relation functions, such as the effective action that governs the dynamics of 
the macroscopic expectation values of the fields. These generating functionals 
contain all relevant physical information about a theory, once the microscopic 
fluctuations have been integrated out. 

The functional RG combines this functional approach with the RG idea 
of treating the fluctuations not all at once but successively from scale to 
scale O [To] . Instead of studying correlation functions after having averaged 
over all fluctuations, only the change of the correlation functions as induced 
by an infinitesimal momentum shell of fluctuations is considered. From a 
structural viewpoint, this allows to transform the functional- integral structure 
of standard field theory formulations into a functional differential structure 
[m [121 [131 [III ■ This goes along not only with a better analytical and nu- 
merical accessibility and stability, but also with a great flexibility of devising 
approximations adapted to a specific physical system. In addition, structural 
investigations of field theories from first principles such as proofs of renormal- 
izability can more elegantly and efficiently be performed with this strategy 

The central tool of the functional RG is given by a flow equation. This flow 
equation describes the evolution of correlation functions or their generating 
functional under the influence of fluctuations. It connects a well-defined initial 
quantity, e.g., the microscopic correlation functions in a perturbative domain, 
in an exact manner with the desired full correlation functions after having 
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integrated out the fluctuations. Hence, solving the flow equation corresponds 
to solving the full theory. 

The complexity of quantum gauge theories and QCD are a serious chal- 
lenge for all field theoretical methods. The construction of flow equations for 
gauge theories has to take special care of the gauge symmetry, i.e., the in- 
variance of the theory under local transformations of the fields in coordinate 
space [m [m [20l EH [22] . The beauty of gauge symmetry is turned into the 
beast of complex dynamical equations and nontrivial symmetry constraints 
both of which have to be satisfied by the flow of the correlation functions. 
Nevertheless, the success of the functional RG in many branches of physics as 
described in other lectures of this volume makes it a promising tool also for 
gauge theories. The recent rapid development of functional methods and their 
application to gauge theories and QCD also in the strong-coupling domain 
confirm this expectation. In combination with and partly complementary to 
other field theoretical methods, the functional RG has the potential to shed 
light on some of the still hardly accessible parameter regions of quantum gauge 
theories. 



2 Functional RG Approach to Quantum Field Theory 
2.1 Basics of QFT 

In quantum field theory (QFT), all physical information is stored in correla- 
tion functions. For instance, consider a collider experiment with two incident 
beams and (n — 2) scattering products. All information about this process 
can be obtained from the n-point function, a correlator of n quantum fields. 
In QFT, we obtain this correlator by definition from the product of n field 
operators at different spacetime points (f{xn) averaged over all possible field 
configurations (quantum fluctuations). 

In Euclidean QFT, the field configurations are weighted with an exponen- 
tial of the action S[ip\, 

{ifixi) . . . ip{Xn)) -.^M j Vififixi) . . . ^{Xn) e~'^^'^\ (1) 

where we fix the normalization A/" by demanding that (1) = 1. We assume that 
Minkowski-valued correlators can be defined from the Euclidean ones by ana- 
lytic continuation. We also assume that a proper regularized definition of the 
measure can be given (for instance, using a spacetime lattice discretization), 
which we formally write as J Vip — > 2?(/?; here, A denotes an ultraviolet 
(UV) cutoff. This regularized measure should also preserve the symmetries of 
the theory: for a symmetry transformation U which acts on the fields, ip — > ip^ , 
and leaves the action invariant, S[tf\ S[lp^] = S[(p], the invariance of the 
measure implies 
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Vip^ Vip^ = / Vip. (2) 

A J A J A 

For simplicity, let (p denote a real scalar field; the following discussion also 
holds for other fields such as fermions with minor modifications. All n-point 
correlators are summarized by the generating functional Z[J], 

Z[J] =c^['^l = / P^e-^[^1+^'^^ (3) 



with source term J J(p = j d^x J{x)p{x). All n-point functions arc obtained 
by functional differentiation: 

{^ij:,)...^i^„))^J-(^^M^) . (4) 
Z[0\ \6J[xi) . . .dJ{xn) J 

Once the generating functional is computed, the theory is solved. 

In Eq. ^ , we have also introduced the generating functional of connected 
correlators^^ W[J] = lnZ[J], which, loosely speaking, is a more efficient way 
to store the physical information. An even more efficient information storage 
is obtained by a Legendre transform of VK[J]: the effective action F: 



r[<^] = sup(^y", 



J<j) - W[J] . (5) 



For any given 0, a special J = Jsup = J[<t'] is singled out for which J Jcj) — 
W[J] approaches its supremum. Note that this definition of F automatically 
guarantees that F is convex. At J = Jsup, we get 



SW[J] 1 SZ[J] 



5 J Z[J] 5 J 



{V)J- (6) 



This implies that (f> corresponds to the expectation value of (p in the presence 
of the source J. The meaning of F becomes clear by studying its derivative 



at J — Jsup 



SF[^]_ fSW[J]SJiy) , 



^Hx) Jy SJ{y) S<j){x) JyS<j){x 

This is the quantum equation of motion by which the effective action F[(j)\ 
governs the dynamics of the field expectation value, taking the effects of all 
quantum ffuctuations into account. 



^ In this short introduction, we use but make no attempt at fully explaining the 
standard QFT nomenclature; for the latter, we refer the reader to any standard 
QFT textbook, such as [231124] . 
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From the definition of tlie generating functional, we can straightforwardly 
derive an equation for the effective action: 

e-^[^l = exp [-S[<l> + ^]+I^J^^y (8) 

Here, we have performed a shift of the integration variable, ip —>■ cp + (f). 
We observe that the effective action is determined by a nonlinear first-order 
functional differential equation, the structure of which is itself a result of a 
functional integral. An exact determination of r[(j)\ and thus an exact solution 
has so far been found only for rare, special cases. 

As a first example of a functional technique, a solution of Eq. ([8]) can be 
attempted by a vertex expansion of F [(/)], 

oo 

m^Y. - / d^^i---d^^n/^<"na;i,...,a;n)0(xi)...<^(x„), (9) 

where the expansion coefficients _r(") correspond to the one-particle irre- 
ducible ( IPI) proper vertices. Inserting Eq. ([9]) into Eq. ([8|) and comparing 
the coefficients of the field monomials results in an infinite tower of coupled 
integro-differential equations for the J^'-"-': the Dyson- Schwinger equations. 
This functional method of constructing approximate solutions to the theory 
via truncated Dyson- Schwinger equations, i.e., via a finite truncation of the 
series Eq. (|9]) has its own merits and advantages; their application to gauge 
theories is well developed; see, e.g., [25l [26j [271 EH]- Here, we proceed by 
amending the RG idea to functional techniques in QFT. 



2.2 RG Flow equation 

A versatile approach to the computation of F is based on RG concepts |14) . 
Whereas a computation via Eq. ((8|) or via Dyson-Schwinger equations corre- 
sponds to integrating-out all fluctuations at once, we can implement Wilson's 
idea of integrating out modes momentum shell by momentum shell. 

In terms of -T, we are looking for an interpolating action 1^, which is also 
called effective average action^ with a momentum-shell parameter k, such that 
/fe for k ^ A corresponds to the bare action to be quantized; the full quantum 
action F should be approached for fc ^ 0, 

Fk->A — 5'baro, A^O = F. (10) 

This can indeed be constructed from the generating functional. For this, let 
us define the IR regulated functional 

e^-W =^fe[J] :=exp(^-A5fe 

VLpe-s'^^]-^s^^v]+! Jv^ (11) 
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where ^ 

A5fcM = i y ^{~q)Rk{qMq) (12) 

is a regulator term which is quadratic in ip and can be viewed as a momentum- 
dependent mass term. The regulator function Rk{q) should satisfy 

lim Rk{q)>0, (13) 

which implements an IR regularization. For instance, if Rk ~ fc^ for ^ fc^, 
the regulator screens the IR modes in a mass-like fashion, ~ k"^. Further- 
more, 

hm Rk{q)^0, (14) 

/q ^0 

which implies that the regulator vanishes for fc — ^ 0. As an immediate conse- 
quence, we automatically recover the standard generating functional as well 
as the full effective action in this limit: Zk^o[J] = Z[J] and /fe_,o — r. The 
third condition is 

lim Rk{q) oo, (15) 

which induces that the functional integral is dominated by the stationary 
point of the action in this limit. This justifies the use of a saddle-point ap- 
proximation which filters out the classical field configuration and the bare 
action, Fk^A 5 -f const.. A sketch of a typical regulator that satisfies these 
three requirements is shown in Fig. [TJ Incidentally, the regulator is frequently 
written as 

Rkip')^p'r{p'/k'), (16) 

where r{y) is a dimensionless regulator shape function with a dimensionless 
momentum argument. The requirements (|13 p -(|15 |) translate in a obvious man- 
ner into corresponding requirements for r{y). 

Since we already know that the interpolating functional Fk exhibits the 
correct limits, let us now study the intermediate trajectory. We start with the 
generating functional Wk [J] , using the abbreviations 

Keeping the source J fixed, i.e., k independent, we obtain 

dtWk[J] = -\j Vipip{-q)dtRk{q)v{q)c-S-^''+f-''^ 

= -\j dtRkiq) Gkiq) + dtASk[<t>]. (18) 

Here, we have defined the connected propagator 

Gkip) = (P) = i'fii-pMp)) - {^{-p)){p{p)). (19) 
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Fig. 1. Sketch of a regulator function Rk{p^) (lower curve) and its derivative 
dtRkip^) (upper curve). Whereas the regulator provides for an IR regularization 
for all modes with p < fe , its derivative implements the Wilsonian idea of integrat- 
ing out fluctuations within a momentum shell near p^ :^ k^. 



(Note that we frequently change from coordinate to momentum space or vice 
versa by Fourier transformation for reasons of convenience.) Now, we are in a 
position to define the interpolating effective action Fk by a slightly modified 
Legendre transform|f| 



AM = sup (^J J<l>- Wk[J]j - ASki 



(20) 



Since we later want to study Fk as a functional of a fc-independent field cj), it is 
clear from Eq. ([20|) that the source J = Jsup = J[(l>] for which the supremum 
is approached is necessarily k dependent. As before, we get at J = Jsup^ 

0(x) = {^{x))j = (21) 
The quantum equation of motion receives a regulator modification. 



From this, we deduc^: 

SJ{x) S^Fki 



Rk{x,y). (23) 



Now, only the "sup" part of 1^ is convex. For finite fc, any non-convexity of F^ 
must be of the form of the last regulator term of Eq. (|20p . 
"* In case of fermionic Grassmann-valued fields, the following derivative should 
act on Eq. (|22p from the right. 
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■dtRk 




Fig. 2. Diagrammatic representation of the flow equation ([28]): the flow of Fk is 
given by a one-loop form, involving the full propagator Gk ~ (r^'^^ + Rk)~^ (double 
line) and an operator insertion in the form of dtRk (flUed box). 



On the other hand, we obtain from Eq. ([2T|) : 

SJ{x') SJ{x')SJ{y) ~ 
This implies the important identity 



Gk{y-x'). 



(24) 



6{x - x') = 



or, in operator notation. 



(5J(a;) 



5^{y) 5J{x') 



d^y(rf [0]+i?fc)(a;,y)Gfe(2/-x'), 



(rf +i?fc)Gfc. 



Here, we have introduced the short-hand notation 



(25) 



(26) 



(27) 



Collecting all ingredients, we can finally derive the flow equation for Fk for 



fixed (j) and at J = Jgup [14] 



dtWk[J]U + / {dtJ)<i) - dt^Sk[<l>] = -dtWk[J]\j - dt^Sk[4>] 



2 

1 



j^^dtRk{q) Gk{q) 



Tr 



(28) 



This flow equation forms the starting point of all our further investigations. 
Hence, let us carefully discuss a few of its properties, as they are apparent 
already at this stage: 

• The flow equation is a functional differential equation for Fk- In contrast 
to Eq. ([U, no functional integral has to be performed to reveal the full 
structure of the equation. 
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We have derived the flow equation from the standard starting point of 
QFT: the generating functional. But a different - if not inverse - perspec- 
tive is also legitimate. We may define QFT based on the flow equation. For 
given suitable initial conditions, for instance, by defining the bare action 
at a high UV cutoff scale k = A, the flow equation defines a trajectory 
to the full quantum theory described by the full effective action F. In 
the case of additional symmetries, the QFT-defining flow equation may be 
supplemented by symmetry constraints to the effective action. 
The purpose of the regulator is actually twofold: by construction, the oc- 
currence of Rk in the denominator of Eq. (|28|) guarantees the IR regular- 
ization by construction. In addition to this and thanks to the conditions 
(fT3|) and (fT4|) . the derivative dtRk occurring in the numerator of Eq. ((28|) 
ensures also UV regularization, since its predominant support lies on a 
smeared momentum shell near ^ k^. A typical shape of the regulator 
and its derivative is depicted in Fig. [T] The peaked structure of dtRk im- 
plements nothing but the Wilsonian idea of integrating over momentum 
shells and implies that the flow is localized in momentum space. 
The solution to the flow equation ([28|) corresponds to an RG trajectory in 
theory space. The latter is a space of all action functionals spanned by all 
possible invariant operators of the field. The two ends of the trajectory are 
given by the initial condition Fk^A — 5'baro, and the full effective action 
A^o = ^■ 

Apart from the conditions (fT3|) - (fT5|) . the regulator can be chosen arbitrar- 
ily. Of course, the precise form of the trajectory depends on the regulator 
Rk- The variation of the trajectory with respect to Rk reflects the RG 
scheme dependence of a non-universal quantity; see Fig. [3l Nevertheless, 
the final point on the trajectory is independent of Rk as is guaranteed by 
Eqs. (fT3|l-(fT5l). 

The flow equation has a one-loop structure, but is nevertheless an exact 
equation, as is signaled by the occurrence of the exact propagator in the 
loop; see Fig. [21 The one-loop structure is a direct consequence of ASk 
being quadratic in the field operator ip which is coupled to the source [29j . 
Perturbation theory can immediately be re-derived from the flow equation. 
For instance, imposing the loop expansion on 1^, = 5 + -|- 

0{h'^), it becomes obvious that, to one-loop order, r'^^'' can be replaced 
by S''-^-' on the right-hand side of Eq. ([28]) . From this, we infer: 



dtRk ( S^^^ + Rk 



-1 



i dt Tr ln(5(2) + i?,) 



The last formula corresponds to the standard one-loop effective action, as 
it should. 
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Fig. 3. Sketch of the RG flow in theory space. Each axis labels a different operator 
which spans the effective action, e.g., 0^, (9(/f>)^, etc. Once the initial conditions in 
terms of the bare action /fc=A ~ -^barc are given, the solution to the flow equation 
l|28p is a trajectory (solid line) in this space of action functionals, ending at the full 
quantum effective action F = Fk^o- A change of the regulator Rk can modify the 
trajectory (dashed lines), but the end point 7^ stays always the same. 



In practice, the versatility of the flow equation is its most important 
strength: beyond perturbation theory, various systematic approximation schemes 
exist which can be summarized under the method of truncations. 

A first example for such an approximation scheme has already been given 
above in Eq. ([5]), the vertex expansion, which now reads 

/"d^a;i...d^a:„r^")(a;i,...,x„)(/.(xi)...0(x„). (29) 

n=0 ' 

Upon inserting this expansion into the flow equation (|28p . we obtain flow 
equations for the vertex functions Z^^"'' which interpolate between the bare 
and the fully dressed vertices. These equations are similar but not identical 
to Dyson-Schwinger equations, as will be discussed in more detail below. 

As a second example, let us introduce the operator expansion which con- 
structs the effective action from operators of increasing mass dimension. Fo- 
cusing in particular on derivative operators, we arrive at the gradient expan- 
sion. For instance, for a theory with one real scalar field, we obtain 



r, = / d^ 



X 



(30) 



where, for instance, Vfc(0) corresponds to the effective potential. 

Further examples of this type can easily be constructed by combining these 
two expansions in various ways. Formally, expansions of the effective action 
should be systematic; this implies that a classification scheme exists which 
classifies all possible building blocks of the effective action and relates them 
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to a definite order in the expansion. Truncations based on such expansions 
should also be consistent in the sense that, once the maximal order of the 
truncated expansion is chosen, all terms up to this order are kept in the flow 
equation. 

Systematics and consistency are, of course, only a necessary condition for 
the construction of a reliable truncation - they are not necessarily sufficient. 
As a word of caution, let us stress that these two conditions do not guarantee 
a rapid convergence of the truncated effective action towards the true result. 
As for any expansion, convergence properties have to be checked separately, 
in order to estimate or even control the truncation errors. 

Only one general recipe for the construction of a truncation exists: let your 
truncation scheme be guided by physics by making sure that the truncation 
includes the most relevant degrees of freedom of a given problem. 

2.3 Euclidean Anharmonic Oscillator 

Let me illustrate the capabilities of the flow equation by a simple example: 
a 0+1 dimensional real scalar field theory, or, in other words, the Euclidean 
quantum mechanical anharmonic oscillator. This system has first been studied 
by RG techniques in [30 . The bare action that we want to quantize is 

S^jdr Qi^ + \u'x' + A,4^ , (31) 

with uj'^, A > 00 We will mainly be interested in the determination of the 
ground state energy which we expect to be predominantly infiuenced by the 
effective potential. Hence, we consider the truncation 

AN = j dT Q.t2 + Vkix)^ . (32) 

For a concrete computation, we have to choose a regulator which conforms to 
the conditions (fT3|) - (fT5|) . The following choice is not only simple and conve- 
nient, it is also an optimal choice for the present problem, since it improves 
the stability properties of our fiow equation 33J , 

i?fc(p) = (fc'-p')^(fc'-p'), (33) 

which implies that dtRk = 2fc^ ^^(fc^ — P^)- On the right-hand side of the flow 
equation we need F^'^^ ~ {—d^ + Vi^'{x))S{T — t'). In order to project the 
flow equation onto the flow of the effective potential, it suffices to consider 
the special case x — const., for which the right-hand side can immediately 
be Fourier transformed, e.g., —d^ — > p^. We obtain the flow of the effective 
potential: 

^ Also the double-well potential with lu^ < Q can be studied with RG techniques, 
see 
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dtVk{x) - ' 



This is a partial difTerential equation for the effective potential. Aiming at 
the ground state energy, it suffices to study a polynomial expansion of the 
potential, 

Vk{x) ^ ^ujIx' + ^Xkx^ + . . . + Ek, (35) 

where the effects of the fluctuations are now encoded in a scale-dependent 
frequency ujk and coupling At; their initial conditions at fc = yl are given by 
the parameters uj and A in the bare action (|3ip . Unfortunately, Ek is not 
identical to the desired ground state energy Eo^k, but differs by further x- 
independent contributions. This becomes already clear by looking at the UV 
limit k A, where the regulator term becomes ~ iyl^z^, contributing a 
harmonic-oscillator-like ground state energy ~ to E. In order to extract 
the true ground state energy from the flow of Ek, 



-Ek = - 
dk TT k'^ -\- Lu 



Ek = - , ,.2 , (36) 



we can perform a controlled subtraction to avoid the build-up of the unphysical 
contributions: in the limit A = cj = 0, the ground state energy has to stay 
zero, -Eo.A; = 0. This fixes the subtraction term for Eq. ((36|) . which then reads: 

d „ 1 



- LO 



Expanding Eq. (f34|) to higher orders, yields 
do 2 fc2 A 



dfc'^'=-"^ {k^ + ulY T' ^-^^^ 

d _24 e /A,y 

where the ellipsis in the last equation denotes contributions from higher-order 
terms ~ x^ , which we neglect here for simplicity. We have boiled the flow equa- 
tion down to a coupled set of first-order ordinary differential equations, which 
can be viewed as the RG (3 functions of the generalized couplings Eo,k,^k, A^. 
These equations can diagrammatically be displayed as in Fig. |4l The diagrams 
look very similar to one-loop perturbative diagrams, but there are important 
differences: all internal lines and vertices correspond to full propagators and 
full vertices (in our simple truncation here, the vertex represents a full run- 
ning Afc and the propagators contain the running w^). Furthermore, one prop- 
agator in each loop carries a regulator insertion, implying the replacement 
G — > GkdtRkGk in comparison with perturbative diagrams. 
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dtEo^k 




Fig. 4. Diagrammatic representation of Eqs. (|37p - (|39p . The diagrams look similar 
to one-loop perturbative diagrams with all internal propagators and vertices being 
fully dressed quantities. One internal line always carries the regulator insertion dtRk 
(filled box). (One further diagram for dt\k involving a 6-point vertex is dropped, as 
in Eq. p9ll .) 

It is instructive, to solve these three equations (P7)) - ([55| with various ap- 
proximations. We begin with dropping the anharmonic coupling A = 0, imply- 
ing that ujk = ^k=A = ^\ integrating the remaining flow of the ground state 
energy yields 

£;oEEi?o.fc=o=^"dfcAiJo , (40) 

which corresponds to the ground state energy of the harmonic oscillator, as 
it should [h = 1). 

Now, let us try to do a bit better with the minimal nontrivial approxima- 
tion: we drop the running of the anharmonic coupling \k — > A, integrate the 
flow of the frequency (l38|) , insert the resulting Wfc into Eq. ([37]) and integrate 
the energy. Expanding the result perturbatively in A, we find 

which can immediately be compared with direct 2nd-order perturbation the- 
ory [31], 

r.PT 1 3 / A \ 105 / A V 

-En = + ^ w + .... (42) 

Our first-order "one-loop" result agrees exactly with perturbation theory, 
whereas the second-order "two-loop" coefficient has the right sign and order of 
magnitude but comes out too small with a ~ 20% error. However, it should be 
kept in mind that we have obtained this two-loop estimate from a cheap cal- 
culation which involved only a one-loop integral with an RG-improved prop- 
agator. 

Of course, we do not have to stop with perturbation theory. We can look 
at the full integrated result based on Eqs. ([57)1 and ([55]) for any value of A. For 
instance, let us boldly study the strong-coupling limit, where the asymptotics 
is known to be of the form. 
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Fig. 5. Left panel: ground state energy of the anharmonic oscillator for w = 2 
versus the anharmonic coupling A: exact result (solid/black line), 2-loop perturbation 
theory (dotted/black line), flow-equation estimate based on Eqs. (|37p . (|38|l (short- 
dashed/blue line) or on Eqs. (|37p - (|39p (long-dashed/red line). Right panel: Relative 
error of the different estimates. 



1/3 , 



(^) ho + o(A-^/^)]. (43) 



The constant ao has been determined in [3S] to a high precision by means of 
large-order variational perturbation theory: ao = 0.66798 . . .. 

With the simplest nontrivial approximation based on Eqs. ((37|) and ((38)) . 



we obtain q,^3^'ES1 _ 0.6920..., which differs from the full solution by 
merely 4%. Now, solving all three equations ([37 |l -([39 |l simultaneously, we find 
^^7j-|39l _ Q gg20 . . ., which corresponds to a 1% error. 

A plot of the ground state energy is depicted in Fig. O (left panel) for lo — 2 
and A = . . . 100. Whereas the perturbative estimate (dotted line) is a good 
approximation for small A, it becomes useless for larger coupling. Already the 
simplest approximation based on Eqs. ([37]) and ([38]) (short-dashed/blue line) 
is a reasonable estimate of the exact result (solid/black line), obtained from 
the exact integration of the Schrodinger equation as quoted in 36J. The flow 
equation estimate based on Eqs. ([57)) - ([5^ (long-dashed/red line) is hardly 
distinguishable from the exact result. For better visibility, the relative errors 
of these estimates are shown in Fig. [5] (right panel). Over the whole range of 
A, the error of our estimate based on Eqs. ([57)l - (|5^ (long-dashed/red line) 
does not exceed 0.3%. 

The quality of the result is remarkable in view of the extremely simple 
approximation of the full flow. In particular, the truncation to a low-order 
polynomial potential does not seem to be justified at large coupling. In fact, 
there is no reason, why the dropped higher-order terms, e.g., ~ , should be 
small compared to the terms kept. 

The lesson to be learned is the following: it does not really matter whether 
the terms dropped are small compared to the terms kept. It only matters 
whether their influence on the terms that belong to the truncation is small or 
large. 
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2.4 Further reading: regulator dependence and optimization* 

The reliability of the solution of a truncated RG flow is an important ques- 
tion that needs to be addressed in detail for each application. In absence of 
an obvious small expansion parameter, the convergence of any systematic and 
consistent expansion of the effective action can be checked by studying the 
quantitative influence of higher-order terms. Whereas such computations can 
become rather extensive, an immediate check can be performed by regulator 
studies. For the exact flow, physical observables evaluated at fc = are, of 
course, regulator independent by construction. However, truncations generi- 
cally induce spurious regulator dependencies, the amount of which provides a 
measure for the importance of higher-order terms outside a given truncation. 
Resulting regulator dependencies of physical observables can thus be used for 
a quantitative error estimate in a rather direct manner. 

Moreover, the freedom to choose the regulator function within the mild 
set of conditions provided by Eqs. p^ - p3|) can also actively be used for 
an optimization of the flow. A truncated flow is optimized if the results for 
physical observables lie as close as possible to the true results; the influence 
of operators outside a given truncation scheme on the physical results is then 
minimized. If a truncated flow is optimized at each order in a systematic 
expansion, the physical results converge most rapidly towards the true result, 
also implying that the optimized flows exhibit enhanced numerical stability. 
Optimization of RG flows has conceptually been advanced in ^37l [38l[33l[39l [8]. 
In particular, an optimization criterion for the IR regulated propagator at 
vanishing field based on stability considerations has lead to the construction 
of optimized regulators for general truncation schemes 33^ ; we have used such 
a regulator already in Eq. (j33p . A full functional approach to optimization 
has been presented in [8]. Loosely speaking, the optimal RG flow within a 
truncation is identified with the most direct, i.e., shortest trajectory in theory 
space. 

Since optimization helps improving quantitative predictions, enhances nu- 
merical stability and can at the same time be used to reduce technical effort, 
it is of high practical relevance. Optimization is therefore recommended for 
all modern applications of the functional RG. 

3 Functional RG for Gauge Theories 
3.1 RG flow equations and symmetries 

Before we embark on the complex machinery of quantum field theories with 
non-abelian gauge symmetries, let us study the interplay between flow equa- 
tions and symmetries from a more general perspective with emphasis on the 
structural aspects. 

Consider a QFT which is invariant under a continuous symmetry trans- 
formation which can be realized linearly on the flelds; let G be the generator 
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of an infinitesimal version of this transformation, i.e., Q(f> is linear in (/). For 
example, a global O(A^) symmetry in a QFT for A^-component scalar fields is 
generated by 

^;"^-/-/d-./(.)^. (44) 



Incidentally, a local symmetry would be generated by the analogue of Eq. 
without the spacetime integral. Together with the invariance of the measure 
under this symmetry, cf. Eq. ([2]), the invariance of the QFT can be stated by 

= Vipge-^+i-'^P. (45) 

In other words, a transformation of the action can be undone by a transfor- 
mation of the measure. 

What does this symmetry of the QFT imply for the effective action? First, 
we observe that Eq. yields 

= lyp^ iy~{gs) + 1 j{g^)^ e-^+z-^^ 

= -(e^)j + e-^M J J(e^)|^^^e^W. (46) 
Performing the Legendre transform, we find at J = Jgup = J[(f\'- 



o = -{gs)j[,] + J ^-^gcb. 



(47) 



The last term is nothing but gr, and we obtain the important identity (Ward 
identity) 

gm^{gs).ji^y (48) 

It demonstrates that the effective action is invariant under a symmetry if the 
bare action as well as the measure are invariant. 

On this level, the statement sounds rather trivial, but it can readily be 
generalized to the effective average action by keeping track of the regulator 
term, 

- {giS + ASk))j[^] - gASk[^]. (49) 

Here, we learn that the whole RG trajectory is invariant under the sym- 
metry if ^S' = and gASk = 0. This is the case if the regulator preserves the 
symmetry. For instance for the globally 0(A^)-symmetric theory, a regulator 
of the form 

= i 1 ^ ^^-p)6^'Rk{p)^\p) (50) 

is such an invariant regulator. 

Note the following crucial point: since the regulator vanishes, ASk 0, 
for fc ^ 0, Eq. (Uni) appears to tell us that gr^^oi't'] = always 
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Fk^A = S + c.t.'s 
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I 

I 
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Fig. 6. Sketch of the RG flow in theory space with symmetries: the symmetry 
relation cuts out a hypersurface W where the action is invariant under the given 
symmetry. A symmetric regulator keeps the invariance exphcitly, such that the RG 
trajectory always stays inside W (solid line). A non-symmetric regulator can still 
be used if non-symmetric counterterms (c.t.'s) are chosen such that they eat up 
the non-symmetric flow contributions (dashed line) by virtue of the Ward identity 
Eq. H51|) . Also in this case, the resulting effective action is invariant, /l~o £ VV. 

holds, implying that the symmetry is always restored at fc = even for a 
non-symmetric regulator. This is indeed true, but requires that the initial 
conditions at fc = yl have to be carefully chosen in a highly non-symmetric 
manner, since 



Here, the last two terms do generally not cancel each other. Therefore, even 
for non-symmetric regulators, the Ward identity tells us how to choose initial 
conditions with non-symmetric UV counterterms, such that the latter are 
exactly eaten up by non-symmetric flow contributions, see Fig. [6l 

In general, it is advisable to use a symmetric regulator, since the space of 
symmetric actions is smaller, implying that fewer couplings have to be studied 
in a truncation. However, if a symmetric regulator is not available, the flow 
equation together with the Ward identity can still be used. In fact, in gauge 
theories there is no simple formalism with a symmetric regulator. 

3.2 Basics of Gauge Theories 

Let us review a few basic elements of quantum gauge theories for reasons of 
completeness and in order to introduce our notation. 
Consider the classical Yang-Mills action, 



GFaIcP] = {GS) + {GASk^A) - GASk=A- 



(51) 




(52) 
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with the field strength 

f;, = d^Al - d^Al + gr'^^A^Al, (53) 

where the gauge field A" carries an internal-symmetry index (color) , and /"^'^ 
are the structure constants of a compact non-abelian Lie group. The hermitean 
generators T° of this group form a Lie algebra and satisfy 

[T", T^] = i/"''=r^ (54) 

e.g., a, 6, c, = 1, 2, . . . , iV^ — 1 for the group SU(iVc). (In the fundamental rep- 
resentation, the T° are hermitean x Nc matrices which can be normalized 
to tr[T°,T*'] = i^°^) 

The naive attempt to define the corresponding quantum gauge theory by 

Z[J] =y"pAe-^™[^l+^^^^2 (55) 

fails and generically leads to ill-defined quantities plagued by infinities. The 
reason is that the measure 2?^^ contains a huge redundancy, since many 
gauge- field configurations are physically equivalent. Namely, the action is 
invariant under the local symmetry 

K^^l- ^M^" + ar'^^'Al ^Al + M^, (56) 

where uj°'(x) is considered to be infinitesimal and differentiable, but otherwise 
arbitrary. The set of all possible transformations forms the corresponding 
gauge group. The generator of this symmetry is 

where = 8^6"^^ + gf'^'"^A'j^ denotes the covariant derivative in adjoint 
representation. It is a simple exercise to show that J d^xoj^Q\A1;^ — (5A°. 
The full symmetry transformation for non- infinitesimal uj°'(x) can be written 
as (A^ ^ A-T-) 

A^^A^^=UA^U-^-^{d^U)U-\ (58) 

with U = U[ui] — e^'^" being an element of the Lie group. Field config- 
urations which are connected by Eq. ([55)1 are gauge equivalent and form the 
gauge orbit: 

[A°;'''']^{A-\A, = Af, C/[H e SU(7Ve)}. (59) 
Here, A'^^^ is a reference gauge field which is representative for the orbit. 
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In order to define the quantum theory, we would hke to dispose of a mea- 
sure which picks one representative gauge-field configuration out of each orbit. 
This is intended by choosing a gauge-fixing condition, 

F"[^]=0. (60) 

For instance, F" ~ c^^^Ji is an example for a Lorentz covariant gauge-fixing 
condition. Ideally, Eq. ([SO]) should be satisfied by only one A° of each orbit. 
(As discussed below, this is actually impossible for standard smooth gauge- 
fixing conditions, owing to topological obstructions [ID]). 

Gauge fixing can be implemented in the generating functional by means 
of the Faddeev-Popov trick which is usually derived from 

1 = JVF^ 5[FT = J dM(c^) S[F^{An] det (^^^) , (61) 

where d/i denotes the invariant Haar measure for an integration over the 
gauge group manifold (at each spacetime point). This rule is reminiscent 
to the corresponding rule for variable substitution in an ordinary integral, 

1 = / df S{f) = J dxS{f{x)) ^ , for /(x) having only one zero. But as 

already this simple comparison shows, Eq. (pT|) is accompanied by the tacit 
assumption that F"^ = picks only one representative and that the Faddeev- 
Popov determinant det(df / 5lo) > 0. Both assumptions are generally not true; 
and both these properties, namely, that several gauge copies on the same 
gauge orbit all satisfy a given standard gauge-fixing condition and that the 
Faddeev-Popov determinant is not positive definite, are characteristic of the 
famous Gribov problem [41 . A direct solution to this problem in the func- 
tional integral formalism is by no means simple; nevertheless, let us assume 
here that a solution exists which renders the gauge-fixed functional integral 
well defined such that we can proceed with deriving the flow equation. We 
will return to this problem later in the discussion of the flow equation. 

As discussed in standard textbooks, it is now possible to show that the 
Faddeev-Popov determinant is gauge invariant, Z\fp[^"] = dct^^-^^-^ — 
Afp[A]. As a consequence, (5[F°[yl]]Z\Fp[^] can be inserted into the fimc- 
tional integral, such that the redundancy introduced by gauge symmetry is 
removed at least for perturbative amplitudes; this renders the perturbative 
amplitudes well defined and S matrix elements are, in fact, independent of 
the gauge-fixing condition. As a result, the Euclidean gauge-fixed generating 
functional, replacing the naive attempt Eq. (|55|) . becomes 

Z[J] = c^W = J VAApp[A]d[f''[A]]e~^™+^'^^. (62) 

The additional terms can be brought into the exponent: 

(5[F"[A]] ^e-^^'^°"^''^"|^_P = e-^«'-[^l, (63) 
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where we have used a GauBian representation of the S functionaL The expo- 
nentiation of the Faddeev-Popov determinant can be done with Grassmann- 
valued anti-commiiting real ghost fields c, c, yielding 

Afp[A] = J VcDce'-l''^"'''''^^'''' = J VcDce'^'^''. (64) 

The ghost fields transform homogeneously, 

c" ^0"+ gr'"'u;''c'', ^ + gf'''"'uj^(f , (65) 
as induced by a corresponding generator, 



e^h^-ffr^M^^TT + ^-^T^)- (66) 



5 c 5 

In perturbative S matrix elements, the ghosts can be shown to cancel unphys- 
ical redundant gauge degrees of freedom. 

Now, we generalize the generating functional by coupling sources also to 
the ghosts, in order to treat them on the same footing as the gauge field, 

Z[J,-q,fi\ = e^I^'"-''! = j VAVcVcc-'^™''^'^^~^'^*-+-f-''^+-f'"'~-l'"'. (67) 

The construction of the effective action F now proceeds in a standard fashion, 

r[A, c, c] = sup (^J JA + J fjc- J C7] - W[J, rj, fy]^ . (68) 

Since F is the result of a gauge-fixed construction, it is not manifestly gauge 
invariant. Gauge invariance is now encoded in a constraint given by the Ward 
identity Eq. (^5]) applied to the present case with the generator 

g^{x) = es(-) + - ^li"^ - gr'^ (c^^ + c-^) . (69) 

Since G'^Sym — 0, the Ward identity boils down to 

W g''r[A,c,c\ - {G'^iS^f + 5gh)) = 0, (70) 

which in the context of nonabelian gauge theories is also commonly referred 
to as Ward-Takahashi identity (WTI). For instance in the Landau gauge, 

F^[A]=d,A'^^, ^-^^ = -d,D';^''[A], (71) 

together with the gauge parameter a ^ 0, cf. Eq. we can work out this 
identity more explicitly by computing ^°5'gf,gh, e.g., in momentum space. This 
is done in Subsect. 13.41 As an example for how the Ward-Takahashi identity 
constrains the effective action, let us consider a gluon mass term, -Tmass — 
\ I ^A^'ii^'ii- It can be shown order by order in perturbation theory that 
the Ward-Takahashi identity enforces this gluon mass to vanish, = 0; for 
more details, see Subsect. 13.41 Hence, the gluon is protected against acquiring 
a mass by perturbative quantum fluctuations because of gauge invariance. 
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3.3 RG flow equation for gauge theories 

From the gauge-fixed generating functional, the RG flow equation for Fk can 
straightforwardly be derived along the lines of Subsect. 12.21 Using a regulator 
term, 

A^fe = -^A;i-p) {Ru,A)tAp)Al{p) 

-c%p){Rt,,^r\p)c\p), (72) 

we obtain the flow equation 

dtFklA, c,c\ = \ TxdtRkAi^k^ + Rk)-^\A - Tr9ii?fe^gh[(/^f ^ + i?fe)"']gh 
= isTrati?fe(rf^ +i?fe)-i. (73) 

The minus sign in front of the ghost term arises because of the anti-commuting 
nature of these Grassmann- valued flelds; the super-trace in the second line 
of Eq. 173]) takes this sign into account. In our notation in Eq. (jTS]) . f'^^ 
is also matrix-valued in field space, i.e., with respect to (A,c, c); therefore, 
[(/^(2) _|_ Jjj,)-!]^ denotes the gluon component of the full inverse of /^^^^ + 
(and not just the inverse of b^r\^l&A8A + Rk,A)- 

Is this a gauge-invariant flow? Manifest gauge invariance is certainly lost, 
because the regulator is not gauge invariant; e.g., at small p, the regulator 
- here being similar to a mass term ~ is forbidden by gauge invariance, as 
discussed abovej^ But manifest gauge invariance is anyway lost, owing to the 
gauge-fixing procedure. Gauge symmetry is encoded in the Ward-Takahashi 
identity. From this viewpoint, the regulator is merely another source of explicit 
gauge-symmetry breaking, giving rise to further terms in the Ward-Takahashi 
identity, the form of which we can directly read off from Eq. (|49| : 

m GRk + GASk - (a(5gf + 5gh + ASk)} = 0. (74) 

Owing to the additional regulator terms, this equation is called modified Ward- 
Takahashi identity (mWTI). With AS'fc being quadratic in the field variables, 
these regulator-dependent terms have a one-loop structure, since (GASk) — 
GASk corresponds to an integral over the connected 2-point function with a 

® One may wonder whether a gauge-invariant flow can be set up with a gauge- 
invariant regularization procedure. In fact, this is an active line of research, and 
various promising formalisms have been developed so far |42II43||44] . However, the 
price to be paid for the resulting simple gauge constraints comes in the form of 
nontrivial Nielsen identities, non-localities or extensive algebraic constructions. 
For practical application, we thus consider the standard formulation described 
here as the most efficient approach so far. 
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regulator insertion ^ R^- Since the standard Ward-Takahashi identity already 
involves loop terms (cf. Subsect. 13. 4p . the solution to Eq. ([73]) is no more 
difficult to find than that of the standard WTI. 
As before, we observe that 

hm Wk = W, (75) 

such that a solution to the mWTI Wk = satisfies the standard WTI, W — 0, 
if the regulator is removed at fc = 0. Such a solution is thus gauge invariant. 
Loosely speaking, the mWTI Wk ~ defines a modified gauge invariance 
that reduces to the physical gauge invariance for k ^ 0. For a discussion 
of these modified symmetry constraints from different perspectives, see e.g., 
[inilinilllllSlESlHlHZlSHlIH], or the review [3]. 

One further important observation is that the flow of the mWTI satisfies 

9^^^ - -\Grmt^G<i- ^^Wk. (76) 

Here, we have used the collective field variable ^ = {A, c, c) . The collective 
indices A, B,C, . . . label these components, and denote all discrete indices 
(color, Lorentz, etc.) as well as momenta; e.g., the flow equation reads in this 
notation: 9* A = ^dtR^^G^"^. The derivation of Eq. ^ from Eq. ^ is 
indeed straightforward and a worthwhile exercise. 

Let us draw an important conclusion based on Eq. (j76[) : if we manage to 
find an effective action Fk which solves the mWTI Wk = at some scale 
k, then also the flow of the mWTI vanishes, dtWk = 0. In other words, the 
mWTI is a fixed point under the RG flow. Now, if this Fk is connected with 
Fk' at another scale k' by the flow equation, also Fk' satisfies the mWTI at 
this new scale, Wk' = 0. Gauge invariance at some scale therefore implies 
gauge invariance at all other scales, if the corresponding I^-'s solve the flow 
equation. The whole concept of the mWTI is sketched and summarized in 
Fig. [3 

Unfortunately, the picture is not as rosy as it seems for a simple practical 
reason: in the general case, we will not be able to solve the flow equation 
exactly. Hence, the identity (|76p and thus dtWk = will be violated on the 
same level of accuracy. This problem is severe if dtWk = is violated by 
RG-relevant operators, see Fig.[8l the latter are forbidden in the perturbative 
gauge-invariant theory. 

In the perturbative domain where naive power-counting holds, an RG rel- 
evant operator is potentially given by the gluon mass term, i J m\A'^A1^. 
Gauge invariance in the form of the standard WTI enforces m\ = as a 
consequence of W = 0, as mentioned above. By contrast, such a bosonic mass 
term in a system without gauge symmetry would receive large contributions 
from fluctuations; perturbative diagrams are typically quadratically divergent 
in such systems. In a naively truncated RG flow of a gauge system, we can 
therefore expect the gluon mass to become large if the gauge symmetry is not 
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/ r^^A = ^-f- c.t.'s 




Fig. 7. Sketch of an RG flow with gauge symmetry in theory space: the standard 
Ward identity W cuts out a hypersurface of gauge-invariant action functionals. A 
gauge-invariant trajectory (solid line) would lie completely within this hypersurface. 
Instead, the presence of the regulator leads to the mWTI Wfe, cutting out a different 
hypersurface, which approaches W in the limit fc — > 0. A solution to the flow equation 
stays within the Wk hypersurface (dashed line), implying gauge invariance of the 
final full action F = /fc=o. 



respected, m\ ^ ■ Now, the mWTI assists to control the situation: since 
gauge symmetry is not manifestly present in our RG flow, wc cannot expect 
the gluon mass to vanish at all values of fc. As discussed in more detail in Sub- 
sect. [Hill the gluon mass becomes of order m\ ~ g^fc^, as can be determined 
from Wfc = 0. As long as perturbative power-counting holds, this implies that 
— > for fc ^ 0, and the gauge constraint becomes satisfied in the limit 
when the regulator is removed. 

This consideration demonstrates that the mWTI can turn a potentially 
dangerous relevant operator, which may appear in some truncation, into an 
irrelevant harmless operator which dies out in the limit of vanishing regulator. 
From another viewpoint, the mWTI tells us precisely the right amount of 
gauge-symmetry breaking that wc have to put at the UV scale A in the form 
of counterterms, such that this explicit breaking is ultimately eaten up by 
the fluctuation-induced breaking terms from the regulator, ending up with a 
perfectly gauge-invariant effective action. 

This simple gluon-mass example teaches an important lesson: for the con- 
struction of a truncated gauge-invariant flow, the flow equation and the mWTI 
should be solved simultaneously within the truncation. As a standard recipe 
[39l [50] , the flow equation can first be used to determine the flow of all inde- 
pendent operators (e.g., transverse propagators and transverse vertex projec- 
tions) which are not constrained by Wfc = 0. Then use the mWTI to compute 
the remaining dependent operators (e.g., gluon mass, longitudinal terms). In 
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rk=A = 5*'+ c.t.'s 




Fig. 8. Sketch of an RG flow with gauge symmetry in theory space: a truncation 
of the effective action introduces an error, implying a possible range of estimates 
for /fc=o as depicted by the extended ellipse. The error can also have a component 
orthogonal to the W or Wfe hypersurface, if the Ward identities are solved on the 
same level of accuracy as the truncated flow. A violation of the mWTI by RG 
relevant operators is particularly dangerous, since it can quickly drive the system 
away from the physical solution (dotted line) and thus must be avoided. 



this manner, the gauge constraint is exphcitly solved on the truncation and 
gauge-invariance of the truncation is guaranteed. It turns out that the use 
of the mWTI (instead of the flow equation itself) for the computation of a 
dependent operator generically corresponds to a resummation of a larger class 
of diagrams [5T] . 

For the remainder of this subsection, let us return to the Gribov problem 
discussed below Eq. ([6T|) . The fact that standard gauge fixings do not uniquely 
pick exactly one representative of each gauge orbit and that the Faddeev- 
Popov determinant hence is not positive makes the nonperturbative definition 
of the functional integral problematic. Any nonperturbative method which is 
related to the functional integral such as the flow equation therefore appears 
to face the same problem. As an example, let us concentrate on the Landau 
gauge and consider the set of all gauge-field configurations that satisfy the 
gauge-fixing condition — d^A'^^ = 0; the Faddeev-Popov operator then is 
SV^IA^]/ 5u)'^ — —dfi_D'^[A\. A perturbative expansion around g ^ goes 
along with the Faddeev-Popov operator at the origin of configuration space, 
—d^D1^[A\ —9^; since this Laplacian is a positive operator, the Gribov 
problem does not play a role in perturbation theory to any finite order. 



An alternative option could be to use only the flow equation together with a 
regulator that does automatically suppress artificial relevant operators. In fact, 
this is conceivable in the framework of optimization [8]. 
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Moving away from the origin, it is useful to consider the following gauge- 
fixing functional, corresponding to the L2 norm of the gauge potential along 
the gauge orbit, 

Ta[uj] = = + 2 / Lo'^d^Al + / Lo^-d^Df)J + 0{lo^), (77) 

J X J X 

where we have expanded of Eq. ([55]) to second order in w. We can identify 
transverse gauge potentials that satisfy the gauge condition as the station- 
ary points of The Gribov problem implies that a gauge orbit does 
not just contain one but many stationary points of We observe that 
the subset of stationary points given by (local) minima of Ta [w] corresponds 
to a positive Faddeev-Popov operator {—d^D^) > 0; this subset constitutes 
the Gribov region fio. Restricting the gauge- field integration to the Gribov 
region, J VA J^^^ VA, cures the most pressing problem of having a poten- 
tially ill-defined generating functional owing to a non-positive Faddeev-Popov 
determinant. Let us list some important properties of the Gribov region, as 
detailed in [52]: (i) I?g contains the origin of configuration space and thus 
all perturbatively relevant field configurations; (ii) the Gribov region is con- 
vex and bounded by the (first) Gribov horizon d^o, consisting of those field 
configurations for which the lowest eigenvalue of the Faddeev-Popov operator 
vanishes; hence Z\fp = on dfio,. 

Does this restriction of the gauge-field integration to the Gribov region 
modify the flow equation? In order to answer this question, let us go back to 
the functional-integral equation for the effective action (without IR regulator) 
in Eq. dH), 

e-^[*l ^ fv^ exp (-S[^ + ^ , (78) 



with the supplementary condition that (1^) = 0, owing to the shifted integra- 
tion variable, cf. Eq. ([5]). Differentiating both sides with respect to (f) yields 



5m ^ / 5S[^ + 4>] 
6(1) \ S(t> 



(79) 

m 



which is a compact representation of the Dyson-Schwinger equations. Note 
that the same equation can be obtained from the following identity: 

0= /p^i_e-^['^+^]+/^^ (80) 
J dip 

which holds, since the integrand is a total derivative. No boundary terms 
appear here, because the action typically goes to infinity, S — > 00, for an 
unconstrained field ip 00. 

Now, the crucial point for a quantum gauge theory with Faddeev-Popov 
gauge fixing is that the identity corresponding to Eq. (I80p holds also for the 
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constrained integration domain J7q. No boundary term arises, simply because 
the Faddeev-Popov operator and thus the integrand vanishes on the bound- 
ary 9J?g- We conclude that the Dyson-Schwinger equations are not modified 
by the restriction to the Gribov region. Finally, the same argument can be 
transfered to the flow-equation formalism by noting that the effective average 
action has to satisfy an identity similar to Eq. ([78]) including the regulator, 

(81) 

with corresponding IR regulated Dyson-Schwinger equations, 

Sirk[^] + ASkm ^ / d{S[ip + (f>] + ASkiv^ + 4^]) 

d(t) \ 5(1) 

The latter can again be obtained from a functional integral over a total deriva- 
tive similar to Eq. (|80p , indicating that a restriction to the Gribov region does 
not modify Eq. ((82|) . The final step of the argument consists in noting that 
the scale derivative dt of Eq. ([5^ yields the flow equation (once differentiated 
with respect to </>). 

To summarize, solving quantum gauge theories by the construction of cor- 
relation functions by means of functional methods (Dyson-Schwinger equa- 
tions or flow equations) precisely corresponds to an approach with a build-in 
restriction to the Gribov region as an attempt to solve the Gribov problem. 
From a flow-equation perspective, the argument can even be turned around: 
taking the viewpoint that the quantum gauge theory is defined by the flow 
equation, we can initiate the flow in the perturbative deep UV where the 
Faddeev-Popov determinant is guaranteed to be positive. Solving the flow, a 
resulting stable trajectory necessarily stays within the Gribov region. 

Let us close this section with the remark that the picture developed so far 
is not yet complete. The restriction to the Gribov region only removes the 
problem of the non-positive Faddeev-Popov determinant. It does not guar- 
antee that we have integrated over the configuration space by picking only 
one representative of each gauge orbit. In fact, even within the Gribov re- 
gion, there are still Gribov copies. Therefore, the integration domain in gauge 
configuration space has to be restricted even further by picking the global min- 
imum of in Eq. (|77p . The resulting space is known as the fundamental 
modular region A. In practice, the explicit construction of A is difficult; for 
instance, finding the global minimum of the gauge-fixing functional !Fa[^] on 
the lattice, corresponds to an extremely involved spin-glass problem. However, 
it has been argued in [52j within a stochastic-quantization approach that the 
problem of Gribov copies within the Gribov region does not affect the correla- 
tion functions and their computation. This stresses even further the potential 
of functional methods for nonperturbative problems in gauge theories. 
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The foUowing subsection is devoted to a detailed discussion of the gauge con- 
straints in the form of the Ward-Takahashi identity (WTI) and its modified 
counterpart in the presence of the regulator (mWTI). Let us start with the 
standard WTI which we derived already in a compact notation in Eq. ([70]) . 



W g'^r[Ac,c] - + 5gh)) - 0, 



(83) 



which represents the gauge-symmetry encoding constraint that the effective 
action F has to satisfy in a gauge-fixed formulation. In order to work out the 
single terms more explicitly, it is useful to go to momentum space; we use the 
Fourier conventions 

A^^ix) = / e'^-A^Kg), c^ix) - / e'«-c'^(g), c'^(x) = / e-^'^^c^q), (84) 

J q J q J q 

where J = J .^^.^ . This imphes for the functional derivatives, for instance, 



-\qx 



^A-Jq) 



etc. 



(85) 



As a result, the generator of gauge transformations t/" reads in momentum 
space 



SA-A-p) 



(86) 



be 



SAUq-p) 



(5c'' (g — p) 



c%q) 



5d'(q — p) 



This allows us to compute the building blocks of the last two terms in Eq. 
in momentum space; we obtain the gauge transforms 

a'^(p)5gf = -p^pMp) - lgf-^-[ Al{q){p - q),{p - q),Al(p - q), (87) 



e"(p)^gh = -gr" / c^{q)p ■ iq+pViq+p) 

Jq 

^ig^ffebfabc I p^^c^q^)A;{p + gi + q2)c^{q2). 
•'91 >92 



Upon insertion into Eq. we arrive at a explicit representation of the WTI 
in terms of full correlation functions in the presence of a source J which is 
field dependent by virtue of the Legendre transform, J = Jsup = J 
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- --gr"' f{p + q)Ap + qUA';X-q)Al{p + g))eo„ 

^ J q 

-gr"' [ pMp + 91 - q2)q2^.5^f - igf^Alip + - 92)) 

x(c'=(gi)c^(92))con 
_.g2fabcffeb f p^{c^{q,)Alip + q,-q2)c/{q2)Un, (89) 

91:92 

where (. . .}con denotes only the connected part of the correlation functions, 
e.g., {(p(p)con = (vv) ~ All terms on the right-hand side are loop terms, 

the last term is even a two-loop term. It should also be stressed that the WTI 
is expressed here in terms of unrenormalized fields and couplings. 

As an example, let us see how the WTI imposes constraints on operators 
in the effective action. For this, we discuss a gluon mass term in the Landau 
gauge, a ^ 0. The gluon-mass operator reads, 

rn..ss^\j y^i\AlAl. (90) 

Its gauge transform yields 

e"(p)r^ass = miip^A^(p). (91) 

This implies that we have to project the remaining terms of the WTI only 
onto the operator ~ p^A'^{p) in order to study the gauge constraint on the 
gluon mass. Let us do so for the first term on the right-hand side of Eq. ([89]) 
by way of example; shifting the momentum q hy q q — p, this term reads 

-gr"' f q^qAAUp-q)At{q))con\c.^o 



a 



-- / Pi^,^,iq)q^GfAp~q,q\A, 

^ J q 



c,c)|a^o, (92) 



where we have introduced the longitudinal projector -Pl.^i^ = qfj.qu/q^ as well 
as the full gluon propagator in the background of all fields G'^i^^j{p — q, q\A, c, c). 
In view of Eq. (|9ip , this propagator is needed only to linear order in the gauge 
field, 

G^,^{JP - q, q\A, c, c) = G^^{p - q, q) 

+Gf,i^{p-q, q-p)V3A,K\p{p-q, -p, q)Gp^{-q, q)Ax{p) 
+C(A2,cc), (93) 

where V^a denotes the full threc-gluon vertex, and all momenta are counted as 
in-flowing. Inserting the order linear in A of Eq. ((93)l into Eq. ((92)) . we observe 
that both gluon propagators are contracted with the longitudinal projector. 
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(94) 



Now, the Landau gauge strictly enforces the gauge fields to be transverse. Any 
longitudinal modes have to decouple in the Landau-gauge limit; in particular, 
we have Gl ~ a ^ 0. As a consequence, the whole expression ([94)) goes 
to zero linearly with a, at least order by order in perturbation theory. We 
conclude that this first term on the right-hand side of the WTI ([89| does 
not support a nonvanishing value for the gluon mass. Let us mention without 
proof that the same property can be shown also for all other terms on the 
right-hand side of Eq. ((89)) . Therefore, the WTI enforces the gluon mass term 
to vanish to any order in perturbation theory, I^mass — 0, implying m\ = 0. 
The WTI protects the zero mass of the gluon against perturbative quantum 
contributions, because of gauge invariance. 

Let us now turn to the modifications of the gauge constraint in the presence 
of a regulator. We already derived the regulator-modified WTI (mWTI) in 
Eq. m. 

Wk + GASk - {GiSgi + 5gh + ASk)) - 0. (95) 

For an explicit representation, we need the gauge transforms of the regulator 
terms, 



g'^{p)ASk,^^ = -5/""= / (f{q)[Rk,Aq + p)-RkM]c\q + p)- 



(96) 
(97) 



Using our previous result for the standard WTI (|89p. the mWTI can now be 
displayed in the more explicit form. 



(^^{d,A^:,f + {d,cnDf{A)c'^ 

-gr"' I {p + q)^{p + qUAl{-q)Al{p + q)) 



1 



-gf 



abc 



P,{Sip + qi - q2)q2^.6'^ ~ '^gf^^J^^ip + qi - 92)) 



x(c'=((7i)c-^(g2))con 
4g2fabcffeb I p^{c^{q,)Al{p + q^~q2)cf{q2))con 

-^gr'' [(i?fc.A)^.(p + q)- {Rk,A)^.iq)] {Ali-q)AUp + q))c 
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gr'^ I [RkAp + q)- RkM] {c%q)c\q+p)) 



con 



(98) 



Jq 



The last two terms denote the modification of the mWTI in comparison to the 
standard WTI. These two terms are one-loop terms with a structure similar 
to the flow equation itself. Both terms vanish in the limit fc ^ and the 
standard WTI is recovered, as it should. Again, we stress that the mWTI is 
expressed in terms of unrenormalized fields and couplings. 

As an example, it would be straightforward to work out the precise con- 
tribution of the regulator terms to the gluon mass which does not vanish in 
contrast to the WTI result [50l|53]. However, here it suffices to estimate the 
order of magnitude of this contribution. The structure Rk(j> + q) — Rkiq), to- 
gether with the fact that we need to project only on the terms linear in A and 
Pf_i (cf. Eq. (|9T|) ). implies that the q integral is peaked around ~ k"^ . The 
dimensionful scales on the right-hand side for the gluon mass are thus set by 
k^ , resulting in in\ ^ g^k^ with a proportionality coefficient that depends on 
the regulator. This is a very unusual bosonic-mass running which guarantees 
that the gluon mass is not a relevant operator in the flow, but vanishes with 



Let us close this subsection by briefly discussing the connection of the 
present formalism with the BRST formalism. The latter involves one further 
conceptual step, emphasizing BRST invariance as a residual invariance of 
the gauge-fixed functional integral. The corresponding symmetry constraints 
on the effective action, the Slavnov- Taylor identities, have the advantage in 
the standard formulation that they are bilinear in derivatives of the effective 
action. This allows for an algebraic resolution of the gauge constraints in 
contrast to the loop computations necessary for the WTI, cf. Eq. (|89|) . If 
the regulator term is included, modified Slavnov- Taylor identities can still be 
derived [HI EH [211 145] , but the result no longer has a bihnear structure. We 
conclude that the BRST formulation has no real advantage in the case of the 
RG flow equation, such that the present formalism can fully be recommended 
also for practical applications. 

3.5 Further reading: Landau-gauge IR propagators* 

The following paragraphs give a short introduction to recent lines of research, 
and may serve as a guide to the literature. 

The functional RG techniques for gauge theories developed above can now 
be used for computing the effective action in a vertex expansion, cf. Eq. (|29|) . 
In momentum space, the expansion reads. 



where = /d^p/(27r)^, S{p) = (27r)^5(-°) (p), and </> = (A;^,c",c"). Insert- 



ing Eq. into the flow equation ([75)1 , we obtain an inflnite set of coupled 



k 



0. 




S{p, + ...+Pn) rl^\pu ■ ■ ■ ,Pn) <P{P1) ■ ■ ■ <P{Pn), (99) 
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first-oder differential equations for the proper vertices /^jj" ■ Truncating the 
expansion at order rimax leaves all equations for the vertices un- 
affected. In order to close this tower of equations, the vertices of order n^ax 
and nmax — 1 can either be derived from their truncated equations or taken 
as bare or constructed by further considerations; see, e.g., [54]. This defines 
a consistent approximation scheme that can in principle be iterated to arbi- 
trarily high orders in rimax. 

Let us consider here the lowest nontrivial order, 

A = iy A1{-q) [ZA{q^)q^ Pxp. + ml5,,]Al{q)^ j c^{q) Z^^M^)q^ ^^(9) + - ■ ■ , 

(100) 

where Pt is the transverse projector, and the ellipsis denotes higher-order ver- 
tices and longitudinal gluonic terms. In the following, we will confine ourselves 
to the Landau gauge a = where longitudinal modes decouple completely; 
moreover, the Landau gauge is known to be a fixed-point of the RG flow 
[49l [46] . The nontrivial ingredients consist in the fully momentum-dependent 
wave function renormalizations Za{p'^) and Zg\i{p^), as well as a gluon mass 
term which has been discussed in detail above. The flow equations for the 
inverse of the transverse gluon and ghost propagators {Gk — {r^^ + Rk)^^), 

4St(p') - Za{p^)p' + ml rl^l^ip') = Zg,(p2)p^ (101) 

are shown in Fig. [9l except for diagrams involving quartic vertices. The flow- 
equation diagrams are reminiscent to those of Dyson- Schwinger equations. But 
there are two differences: Dyson- Schwinger equations also involve two-loop 
diagrams, whereas the flow equation imposes its one-loop structure also on the 
propagator and vertex equations. Second, all internal propagators and vertices 
are fully dressed quantities in the flow equation, whereas Dyson-Schwinger 
equations always involve one bare vertex. Both vertex expansions are slightly 
different infinite-tower expansions of the same generating functional, with the 
flow equations being amended by the differential RG structure. 

A truncation at this order requires information about the triple and quar- 
tic vertices. In a minimalistic approach, they may be taken as bare (possibly 
accompanied by fc-dependent renormalization constants). In general, this non- 
perturbative approximation is expected to be reliable at weak coupling. For 
instance, the perturbative result is rediscovered at high scales, k = Ap] for 
momenta p^ larger than this perturbative scale A^, the wave function renor- 
malizations yield, 

Za{p)^Za^,a {^-^A^^i^-^J. 

Z^y^ip^) Za^.sI. (l - %h ^^^^ In J) , (102) 
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Fig. 9. Flow equations for gluon and gliost propagators in a vertex expansion. 
All internal lines and vertices denote fully dressed quantities, indicated by filled 
circles. To each diagram, there is another one with identical topology but with 
the regulator insertion occurring at the opposite internal line. The ellipses denote 
diagrams involving quartic vertices which are not displayed. 

where r]A ~ —13/22 and rygh — —9/44 denote the anomalous dimensions for 
gluons and ghosts, respectively. Za^.a, ZA^^gh are the normalizations of the 
fields, and aAp is the value of the coupling constant at yip. 

At first glance, there is no reason why the higher-order vertex structures 
which are dropped in the minimalistic truncation should not become domi- 
nant at strong coupling. However, as we learned from the anharmonic oscil- 
lator example in Subsect. 12.31 the quality of a low-order truncation does not 
depend on how large higher-order terms may get but whether they exert a 
strong influence on the low-order equations. Moreover, mechanisms may exist 
that systematically suppress higher-order contributions, such as an IR sup- 
pression of certain propagators; since higher-order vertex equations involve 
more propagators, such a propagator suppression would control large classes 
of diagrams. 

In recent years, evidence has been provided that such a suppression is 
indeed operative in the gluon sector in the Landau gauge: low-order vertex 
expansions reveal an IR suppressed gluon propagator which renders the con- 
tributions from higher gluonic vertices subdominant. This solution has been 
pioneered in [SI] using truncated Dyson-Schwinger equations, see [151 US] for 
reviews; IR gluon suppression in the Landau gauge has meanwhile been con- 
firmed by many lattice simulations [51 [571 UHl UHl [SnUSl [S1IS3]- In the 
continuum, this scenario goes along with an IR enhanced ghost propagator, 
i.e., IR ghost dominance. This IR enhancement does not spoil the vertex ex- 
pansion owing to a nonrenormalization theorem for the ghost-gluon vertex 
|64| . the running of which is thus protected against strong renormalization 
effects. Also ghost dominance has been observed on the lattice (651 [61], even 
though Gribov-copy and/or finite- volume/size effects appear to affect the IR 
ghost sector more strongly. 
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The nonrenormalization theorem of the ghost-gluon vertex in the Landau 
gauge gives rise to a nonperturbative definition of the running coupling in 
terms of the wave function renormalizations, 

In the IR, the above-described scenario of ghost dominance and gluon sup- 
pression is quantitatively observed in terms of a power-law behavior of the 
wave function renormalizations H 

za{p') - ^gh(p') - ip'r, (104) 

where k denotes a positive IR exponent. In the functional RG framework, 
this solution can be shown to be a fixed point of the flow equations for the 
propagators, cf. Fig. [9l in the momentum regime <^ ^qcd ES]- The 

interrelation of the ghost and gluon propagators owing to the simultaneous 
occurrence of the exponent k arises in all functional approaches from self- 
consistency arguments; as a direct consequence, the running coupling (I103P 
approaches a fixed point in the IR, a(p^ <C ^qcd) ~^ '^*- For instance, a 
truncation with a bare ghost-gluon vertex results in k ~ 0.595; possibly in- 
duced momentum dependencies of the vertex can lead to slightly lower values 
0.5 < K < 0.595 [B7]. Regulator dependencies arising in an RG calculation 
also lie in this range [55] • The IR solutions (|104p are IR attractive fixed-point 
solutions for a wide class of initial conditions and momentum-dependencies of 
the gluon vertices, once the gluon propagator has developed a mass-like struc- 
ture at intermediate momenta at a few times ^qcd [53] . Using suitable vertex 
ansatze, full solutions connecting the perturbative UV branch Eq. (jl02l) and 
the IR power-laws (I104p have been found with Dyson-Schwinger equations, 
see ISH mi- 
Most importantly, the IR power-law behavior featuring gluon suppression 
and ghost dominance agrees with criteria which are expected to be satis- 
fied in two different scenarios of confinement: the Kugo-Ojima [70] and the 
Gribov-Zwanziger [411 1711 [7^ confinement scenario. In particular, a strongly 
IR divergent ghost propagator represents a signature of confinement in these 
scenarios, which describe the absence of color charged asymptotic states and 
(indirectly) a linear rise of the potential between a static quark-antiquark 
pair. The study of correlation functions in connection with these scenarios 
of low-energy gauge theories clearly demonstrates the potential of functional 
methods to access even the strongly-coupled gauge sector by analytical means. 



In the Dyson-Schwinger literature, the gluon and ghost propagator behavior is 
often characterized by dressing functions ZdsEjGdse which are related to the 
wave function renormalizations by Za{p^) ~ ■Z'dse(p^) and Zgh(p^) — G^\^{p'^) 
for fc -> 0. 
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4 Background-field flows 

Imagine for a second that we knew nothing about computing loops and con- 
structing amphtudes in some sort of expansion which involves a perturbative 
or even a fully dressed propagator. If we knew only the degrees of freedom 
of our gauge system and the symmetries we would be trying to write down 
an effective action in terms of all possible gauge-invariant gluon operators, 
such as F^^F^^,, F^jj{D ,iY^ F^^ etc. and determine the coefficients, e.g., in 
a manner similar to chiral perturbation theory. The result would be gauge in- 
variant by construction, and we would never worry about Ward identities and 
how gauge invariance can be encoded in a nontrivial manner in a gauge-fixed 
formulation. 

4.1 Background-field formalism 

The background-field formalism aims precisely at the construction of such an 
effective action, nevertheless by computing loops and integrating out fluctua- 
tions in a special gauge-fixed manner. Here is a rough sketch of the idea: 

(1) Introduce an auxiliary, non-dynamical field (background field) with its 
own auxiliary symmetry transformation Q. 

(2) Construct a gauge-fixed QFT which has broken invariance under the stan- 
dard gauge transformation Q but a manifest invariance under {Q + Q). 

(3) Let the full F inherit the symmetry properties in the end by setting A — A 
after the gauge-fixed calculation. 

Let us start with (1): we introduce A° and corresponding covariant derivatives 
7?°^ = dfiS""^ + gf°''*A'jj^, and the generator of symmetry transformations 



which we call the background transformation. Note that the ghosts are not 
affected by Q. Together with Eq. ([69|) . it is obvious that {A— A) now transforms 
homogeneously under Q + Q, 



As step (2), we choose a gauge fixing which fixes the Q symmetry but is 
invariant under G + Q: 




(105) 



/ 



d^ycj\y)ig + g)\y)iA 



Ar^ix)^gr'^u;\x)iA-A)lix). (106) 




(107) 



In fact, the gauge- fixing term is invariant under the combined transformation. 
With the Faddeev-Popov operator 
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Xpa. 

= -Dl^D^^, (108) 

it is also straightforward to show that 

{Q + ^)5'gh = -{Q + Q) j d^:r ?fDl^D''^c'' = 0. (109) 

Finally, let us merely sketch step (3): the price to be paid so far is that F 
now depends on A and A. But at the end of the calculation, we can identify 
A = A, such that 

o^{g + g)r[A,A]\^^^^gr[A,A], (no) 

where the first equality holds by construction and the second arises from 
setting A = A. Now, it is possible to prove that the background effective action 
with A = A reduces precisely to the standard effective action, r[A, A] = 
r[A] [73l EH in]. Therefore, Eq. (jUOp verifies the desired gauge-invariant 
construction of r[A]. This r[A] thus only consists of gauge-invariant building 
blocks. Of course, there is still a nontrivial constraint which becomes visible 
if we go away from the limit A — A; namely, Qr[A,A] has to satisfy the 
standard WTI 07]. 



4.2 Background-field flow equation 

The desired properties of the effective action expressed by Eq. (I110|) can be 
maintained in the construction of the corresponding RG flow equation, if also 
the regulator satisfies 

(^ -t- ^)A5fc = 0. (Ill) 

This holds, e.g., for the choice 

ASk = lji^- A)RkA^A){A j ci?fe,gh(4h)c, (112) 

where Z\A.gh are operators that can depend on A and transform homoge- 
neously. For the gluon sector, a suitable choice can, for instance, be given 

by 

i^A^l = {-DfD':S,., -DfDi^S,, + 2ig{Fl,T'r, ■ ■ ■}, (US) 

where the first form corresponds to the background-covariant Laplacian, and 
the second also contains the spin-one coupling to the background field. For the 
ghost sector, the Laplacian is also an option, = —D^^'D'^^^. The resulting 
flow equation in the background-fleld gauge reads [75] 

dta[A,c,c,A] = isTr {dtRk{A)[rl^^\A,A]+ RkiA)]-^Y (114) 
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Here, it is a temptation to set A — A in the flow equation; however, the above 
construction tells us that this should be done only at the end of the calculation 
at A; = 0. 

Let us parameterize (suppressing ghosts for a moment) the effective action 
as [75] 

A] = rr[A] + rf A], (115) 

where -T™^ is a gauge-invariant functional, and r^^'^^° denotes the gauge- 
non-invariant remainder that satisfies r'l''"*^'' [A, A = A] = 0, cf. Eq. pTU)) . 
Considering the second functional derivative with respect to A, 

rl^'^ [A, A] = r;""(') [A] + rf [A, A], (ii6) 

we observe that must be singular. This is because of the zero modes 

associated with gauge invariance: a variation with respect to A which points 
tangentially to the gauge orbit has to leave F™"^ invariant, corresponding to a 
flat direction. For the flow equation pi4p to be well defined, the contribution 
of /^s™se(2) ^j^g denominator in Eq. (I114p has to lift these flat directions. 
This demonstrates that /^^augc j^^g^- dropped from a truncation, even 

though we may ultimately be interested only in F™"^ . 

On the other hand, the mWTl does not impose any constraint on ^ 

= Wfc[A] = Wfc[rr"1: (117) 

which implies that we have the full freedom to choose any gauge-invariant 
functional as an ansatz for F™^ , and its solution will solely be determined by 
the flow equation. 

To summarize, the background formalism facilitates the construction of a 
gauge-invariant RG flow in which the manifestly gauge-invariant parts of the 
effective action r'™^ can be separated from the gauge-dependent parts r^'^^^'^ . 
In practice, we can construct our ansatz for F™^ by picking gauge-invariant 
building blocks, in a similar manner as we would do for other effective field 
theories. In addition, we have to construct (a truncation for) psauge .^j^-j^ -j-j^g 
aid of the mWTI ()117|) . in order to lift the gauge zero modes in the flow 
equation. 

4.3 Running coupling 

The background formalism also provides for a convenient nonperturbative 
definition of the running coupling. As a general remark, let us stress that 
there is no unique definition of the running coupling in the nonperturbative 
domain. Universality of the running coupling holds only near fixed points; e.g., 
only the one-loop coefficient of the perturbative (3 function is definition and 
scheme independent (in a mass-independent regularization scheme, also the 
two-loop coefficient is universal). Hence, any result for the running coupling 
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in the nonperturbative domain has to be understood strictly in the context of 
its definition. 

The definition within the background formalism follows, for instance, from 
the RG invariance of the background-covariant derivative, which is a gauge- 
covariant building block of F™^, 

D;''[A]=d^5^' + gr''^A';,. (118) 

Here, we have used the notation g for the bare coupling. Obviously, the first 
term 9^5"* is RG invariant. Hence, also the product of g and A must be RG 
invariant 173 ] , 

dt{gAl)^0, ^ gA^^gA^^^^, (119) 

where g now denotes the renormalized coupling and A^ the renormalized 
background field. Renormalization of the background field is described by a 

~ 1/2 ~ 

wave function renormalization factor, A^ = A. Consequently, also the 
running of the coupling is tied to the same wave function renormalization, 

= z^'f- (120) 

We obtain for the (3 function of the running coupling 

/3g2 = dtg"" = w^ 1] - -dt InZfc, (121) 

where rj denotes the anomalous dimension of the background field. The wave 
function renormalization of the background field can be read off from the 
kinetic term of the background gauge potential, 

rr[A]^ I + (122) 

where the dots represent further terms in the truncation. The running coupling 
is thus linked to the lowest-order term of an operator expansion of the effective 
action. 

According to its definition, the running coupling can be viewed as the re- 
sponse coefficient to excitations about the background field. If the background 
field, as a natural choice, is associated with the vacuum state the running cou- 
pling measures the coupling between the vacuum and excitations, for instance, 
in the form of (effective) quark and gluon fluctuations. 

4.4 Truncated background flows 

Let us study the background flow in a simple approximation, following '76j. 
In particular, we will be satisfled by a minimal truncation for J^^augc^ p^^. ^j^jg^ 
we expand I^saugc leading order in A ~ A, 



rr^''[A,A] = j{A^ A)l m;1 {A - At + 0{{A - Af). (123) 
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Then, we fix M^^ by the tree-level order of the mWTI: 



jngaugo 



A,A] = — liA- A)l i-Dl^Df) {A - A)l 



(124) 



which is just the classical gauge-fixing term. This approximation has an im- 
portant consequence: to this order, 



r, 



(2) 



r 



42) 



0(2) 



(125) 



k — ^ k "T i J, 

is independent of {A — A)^ and we can set A — A under the flow for all k. The 

- fe 



form of r^^^^° in Eq. (|124p is just enough to lift the gauge zero modes. This 



gives us an approximate flow for Ff 
dtrriA] = I Tr I dtRkA^^ 



inv 

k ' 



r, 



inv(2) 



-i-DD) + Rk{A) 



-Tr IdtRkM^sh) 



r, 



inv(2) 



Rk{A) 



gll 



(126) 



where we have dropped the bars on the right-hand side, since A = A. So 
far, we have suppressed the ghost fields c, c. If the dependence of r'™'^ on the 
ghosts is such that ghosts are contracted with homogeneously transforming 
color tensors, the mWTI is satisfied to the same level of accuracy as by the 
ghost-independent part. The classical ghost action in background-field gauge 
reads 5gh = — / xc"^D°'^D^'^c^\ hence, the lowest-order approximation for 
our truncation is given by the classical term at A — A, 



7-iinv 



(127) 



Let us now concentrate on the running coupling in the background gauge. 
As discussed above, this can be read off from the kinetic terms of the gauge 
field which we will therefore choose as the only nontrivial part of our simplest 
nontrivial truncation, 



T^inv 
^ fe.gli 



[A] 



rpa rpa 



(128) 



where the running of the wave function renormalization remains to be calcu- 
lated. This is done by projecting the right-hand side of the fiow equation onto 
this kinetic operator only; any other operator generated by the flow, e.g., con- 
taining ghost fields, is dropped. Hence, we can set c, c = in F^^'' (of course, 
after functional differentiation), implying that I^^^-* becomes block-diagonal 
with respect to gluon and ghost sectors. In the gluon sector, we obtain 



^xSA; 



cb 



SAl (129) 
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where we have introduced the notation 



(130) 



for the covariant spin-one Laplacian. Using the Feynman gauge a = 1 here 
for simplicity, the gluon sector thus reduces to 



/ pinv(2)xafc | _ 7 7-)^- 



,ab 



(131) 



From Eq. ()127|) . we can immediately read off the form of F^^'l in the ghost 
sectorial 



inv(2)Nafc| 



-\cb 



Upon insertion into Eq. (|126p . the flow equation boils down to 



(132) 



dtrriA] - - Tr [dtRkA^kT^T + Rk,Ay'] -Tr [dtRk,sh{-D^ + Rk,sh)~^] ■ 

(133) 



A convenient choice for the regulator is given by (cf. Eq. ([16 



Rk,A = Zk Vt r(I?T/fc'), i?fc,gh = i-D') r{-D^/k^), lim r(y) ^ -. 

y^O y 

(134) 

The insertion of the wave function renormalization Zk into the regulator is 
useful, since it maintains the invariance of the flow equation under RG rescal- 
ings of the fields (the ghost wave function renormalization has been set to 
Zk,g)x = 1 in our truncation). The choice r{y) — > ^ implies that the IR modes 
with < k"^ are regulated by acquiring a mass term ~ fc^; the identification 
of the one-loop running is more straightforward with this choice. Then, the 
flow equation reads 



dtrr[A] = \ tt 



dt{ZkriVT/k^)) 



- Tr 



dtri-D'^/P) 



Zfe(l + r(2?T/fc2))J ^' [l + r{~D^/k^) 
=: TvHz^VT/k^) - 2TTni-D^/k^ 



(135) 



The obvious definition of the 7i functions in the last line expresses the fact 
that both terms on the right-hand side are traces over a function of a single 
operator. It is useful to formally introduce the Laplace transform of the Ti, 
functions, 

poo 

(136) 



Jo 

such that the flow equation can be written as 



^ Be aware of footnote 2] on page[71 
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/•CXD /'OO 

Jo Jo 

(137) 

As a result, we have brought the flow equation into propertime form. Note 
that this was possible because of the particular operator structure resulting 
from our truncation: we were able to choose the regulators such that they 
depend on the operators which appear as entries of the block-diagonal F^^-* . 
Such a propertime form of the flow can, for instance, always be established 
to leading order in a derivative expansion of the RG flow [77] . In fact, a wide 
class of truncations can be mapped onto a propertime form [29l [78] which is 
computationally advantageous; moreover, propertime flows have extensively 
and successfully been used in the literature [79 t [80 l [8T 1 [82 t [83 l [84], In the 
present case, we are finally dealing with traces over operator exponentials, so- 
called heat kernels, which can conveniently be dealt with by standard methods. 
The heat kernels of the covariant Laplacians occurring above have frequently 
been studied in the literature; see, e.g., [78]; here we merely need the term 
quadratic in the field strength. 



(138) 



3 (47r)^/2 \kV J 4 A"" 



Let us from now on concentrate on four dimensional spacetime, D = A. From 
Eq. ()137|) . we can extract the flow of the wave function rcnormalization 

+ — (4^X 

- 3 (47r)2^^^^^) + — (i^^^'^)' 

where we have used Eq. (I136|l for y = 0. Together with dt{Zkr{'DT /k"^)) = 
— Zk[2'DT / k'^r' (Vt / k"^) + ?/r(PT/^^)], and the anomalous dimension rj = 
—dt\nZk as defined in Eq. p2ip . we observe that the H functions for zero 
argument uniquely yield 

Wz.(0)-1-|, H(0) = 1, (141) 

independently of the precise form of the regulator shape function r{y) intro- 
duced in Eq. (|134p . This is a direct manifestation of the regularization-scheme 
independence of the one- loop /3 function coefficient, as will become clear soon. 

Introducing the renormalized coupling = Z/^^g"^, cf. Eq. (|120p . we ob- 
tain 

(142) 

' 3 (47r)2 3 (471)2'^- ^ ^ 
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This brings us to our final result for the Yang-Mills [3 function 

Pg2=dt9 =W = ^7^(1 —JJZ^] (143) 



3 (47r)2 V 3 {Any 
22Nc 9^ 220N^ 
3~ (47r)2 9 (47r)4 



(144) 



It is instructive to compare our result to the full perturbative two-loop /3 
function: 

" ~^Ja^ 9 (4^)4 ^^^^^ 

Obviously, we have rediscovered the one-loop result exactly, as we should. 
Moreover, the two-loop coefficient conies out remarkably well with an error of 
only 8%. This is astonishing in two ways: first, we have dropped a number of 
operators that contribute to the two-loop coefficient coupling, e.g., {F^u^ilvY 
{F^F^Y' : ^'^d also the ghost wave function renormalization, which hence 
appear to be less relevant. Second, whereas Eq. ()145p is a universal result 
obtained within a mass- independent regularization scheme such as MS, our 
result arises from a mass-dependent regularization scheme with a mass scale fc; 
for such a scheme, the two-loop coefficient generally is not universal. However, 
in the present truncation, our [3 function indeed is universal to all orders 
computed above, owing to the fact that the dependence of the regulator shape 
function drops out in Eq. (|14ip . We conclude that this simple truncation 
contains already much relevant information about the universal part of this 
two-loop coefficient. Larger truncations indeed contribute further terms which 
are partly non-universal, as expected. An exact two-loop calculation based on 
the functional RG can be found in [85j . 

In view of the quality of this simple truncation, it is tempting to specu- 
late whether the result also contains reliable nonperturbative information. 
However, our result for the (3 function in Eq. (|143p develops a pole at 
g^ = 3(47r)2/(10A'c), clearly signaling the breakdown of the truncation. This 
could already have been expected from the physics content of the trunca- 
tion: independently of the coupling strength, the resulting quantum equations 
of motion are identical to classical Yang-Mills theory. This equations allow 
for plain-wave solutions, describing freely propagating gluons. The trunca- 
tion misses the important IR phenomenon of confinement and thus cannot 
be reliable in the deeply nonperturbative domain. For a stable flow from the 
UV down to the IR, a larger truncation is required that in addition has the 
potential to describe the relevant degrees of freedom for confinement. 



4.5 Further reading: IR running coupling* 



The nonperturbative estimate of the (3 function for the running coupling de- 
rived above has a remarkable property. All higher-loop order corrections can 
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be summed into a geometric series, resulting in the structure of Eq. (jl43p . The 
origin of this structure hes in the fact that we have included the wave function 
renormalization in the regulator, cf. Eq. (I134p . The apparent formal reason 
was that the RG invariance of the flow equations is thus maintained; but in 
addition, we thereby obtain a result which contains a resummation of a larger 
class of perturbative diagrams. 

From a more physical viewpoint, the wave function renormalization de- 
scribes the deformation of the perturbative spectrum of fluctuations, 5^^^ ~ 
, as it is induced by quantum fluctuations; at lower scales, we find the spec- 
trum r'j^^ ~ ZkP^. Therefore, the inclusion of in the regulator leads to a 
better adjustment of the regulator to the deformation of the spectrum, i.e., 
to the spectral flow of Ff, . 

If Zk had not been included in the regulator, we would have found only the 
one-loop (3 function in the truncation (|128p without any higher-loop orders. 
The latter would have been encoded in the flow of higher-order operators 
outside this simple truncation. For instance, the inclusion of the higher-order 
operator {F^yF'^^Y would have resulted in a /3-function estimate of the form 
/3g2 = -qg"^, with 

with some coefficient &i, and being the correct one- loop result, and W2 
denoting the generalized coupling of this higher-order operator. In this trun- 
cation, all nonperturbative information is contained in the flow of W2 , which in 
turn can reliably be computed only by including even higher-order operators. 
A good estimate therefore probably requires a very large truncation. Even if 
the precise infrared values of the higher couplings u;2, W3, . . . may not be very 
important, their flow exerts a strong influence on the running coupling in this 
approximation. 

Together with the inclusion of Zk in the regulator, our estimate for the /3 
function would be of the form (3g2 = ■qcf' with 



'(4^ ' "-^ (47r 



(147) 



with a further coefficient ^2 , and di < can be read off from Eq. (|143p . Par- 
ticularly this di makes an important contribution to the two-loop (3 function 
coefficient, as discussed above. Contrary to Eq. (|146p . this equation contains 
information to all orders in g^, even for the strict truncation W2 = 0, solely 
due to the spectral adjustment of the regulator. 

Of course, for more complicated truncations, the deformation of the spec- 
(2) 

trum of F^ becomes much more involved. A natural generalization would 

(2) 

thus be the inclusion of the full F^, ' in the regulator. However, any depen- 
dence of Rk on the fluctuation would invalidate the derivation of the flow 
equation in Eq. (|18p. especially spoiling the one- loop structure. But within 
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— (2) 

the background formalism, we can at least include -T^ where we have set all 
field dependence of F equal to the background field A, 

) = ^f K^f /(^fcfc'), (148) 

such that the regulator is fully adjusted to the spectral flow of the fluctuation 
operator evaluated at the background field. The resulting flow equation then 

— (2) ~ (2) ~(2) 

contains also P^; derivatives; for instance, a,t A = A where rj: ' = , we 
find 



dtFk = - STr 
2 



(2) 

-yr'iy) , dtF^ r{y)+yr'{y) 



1 + r{y) r'^) 1 + r{y) 



(149) 



Despite the additional terms, this form of the fiow equation (together with 
the approximation of setting A = A already for finite k) is technically advan- 
tageous, since it can be written in generalized propertime form by means of 
a Laplace transformation, see the discussion in the preceding section. 

(2) 

The additional dtFf. terms are the generalization of the rj term in 
Eq. ()141|1 and thus can support a further resummation of a large class of dia- 
grams. In this has been used together with an operator-expansion trun- 
cation involving arbitrary powers of the field strength invariant {F^^F^^)'^ , 
n — 1,2,..., to get a nonperturbative estimate of the (3 function. The anoma- 
lous dimension has the structure. 



l + di 



^ - , (150) 



(47r)2 ^ {AtiY ^ (47r) 



the form of which can approximately be resummed in a closed-form inte- 
gral expression, see [78]. The resulting (3 function exhibits a second zero at 

^ gl > 0, corresponding to an IR fixed point, see Fig. [TOl Hence the 
flow-equation results for the IR running coupling in background gauge based 
on an operator expansion show strong similarities to those in the Landau 
gauge based on a vertex expansion mentioned in Subsect. 13.51 Within low- 
order perturbation theory, the universality of the running coupling is a well- 
known property. Independently of the different definitions of the coupling, 
the one-loop (and in mass-independent schemes also the two-loop) /3-function 
coefficient is always the same. Here, we also observe a qualitative agreement 
between the Landau gauge and the background gauge in the nonperturbative 
IR in the form of an attractive fixed point. This points to a deeper connection 
between the two gauges which deserves further study and may be traced back 
to certain non-renormalization properties in the two gauges. 

The background formalism has also been used for finite-temperature stud- 
ies of Yang-Mills theory and the approach to chiral symmetry breaking, 
|86[ 187] , and for a study of nonperturbative renormalizability of gauge theo- 
ries with extra dimensions [88^. The background- field formalism lies also at 
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Fig. 10. Running coupling Qs versus momentum scale k in GeV for gauge group 
SU(2), using the initial value as{Mz) — 0.117 for illustration. The solid line repre- 
sents the result of an infinite-order resummation of Eq. (|150p as taken from ;78_ in 
comparison with one-loop perturbation theory (dashed line). 

the heart of a series of RG flow-equation studies of quantum gravity 89J; it 
is, of course, also useful for the study of abelian gauge theories such as the 
abelian Higgs model [301 or strong-coupling QED [21] . An alternative strategy 
to exploit a background field in RG flow equations has been proposed in [52] • 
Finally, the spectral adjustment of the regulator described here is a special 
case of the general strategy of functional optimization [5], cf. Subsect. 12.41 

5 From Microscopic to Macroscopic Degrees of Freedom 

A typical feature of strongly interacting field theories is given by the fact that 
macroscopic degrees of freedom can be very different from microscopic degrees 
of freedom. For instance in QCD, quarks and gluons represent the microscopic 
degrees of freedom, whereas macroscopic degrees of freedom are mesons and 
baryons. The latter are bound states of quarks and gluons. Prominent repre- 
sentatives are the light pseudo-scalar mesons (pions, kaons, . . . ) which carry 
bi-fermionic quantum numbers, ~ ^pip. This type of fermionic pairing occurs 
in many systems, also in condensed-matter physics with strongly correlated 
electrons. Generically, a strong fermionic (self-)interaction is required for this 
pairing. In QCD, this is, of course, induced via the interactions with glu- 
ons. For an efficient description of the physics, it is advisable to take this 
transition from microscopic to macroscopic degrees of freedom into account 
[93l[Mlll3IMl9Zll98]. 
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5.1 Partial bosonization 

An explicit example for a fermion-to-boson transition is provided by the 
Hubbard-Stratonovich transformation, or partial bosonization. Let us discuss 
this transformation with the aid of a specific system: the Nambu-Jona-Lasinio 
(NJL) model [99] . We consider a version with one Dirac fermion, defined by 
the action 

5f = / d^^: (^i# + I Xii^Pi:)^ - (V^Ts^)'] V (151) 



This model has a U(l) xU(l) symmetry which here plays a similar role as chiral 
symmetry in QCD; in particular, it protects the fermions against acquiring a 
mass due to fluctuations. 

Partial bosonization is obtained with the aid of the following mixed 
fermionic-bosonic theory, 

Sfb = J d*x {i/ii^V + + /^[V^^'l^V' - ^^'r0>]} , (152) 

with Pl.r = (1/2)(1 ± 75) being the projectors onto left- and right-handed 
components of the Dirac fermion. In fact, the model (|152p is equivalent to 
that of (|15ip also on the quantum level if 

For a proof, it suffices to realize that 

^j^^-l{^i9'4'+\7^[(^^)'-(i'^^^^f]) ^AAe-^^ (154) 

where we have used that {■ipPj^tp)(-ipPF(,tp) = i'ipip)'^ — ('075^')^; and JV abbrevi- 
ates the Gaufiian integral over which is a pure number and can be absorbed 
into the normalization of the remaining fcrmionic integral J 'Dip'D^p e~^^ . Also 
on the classical level, the equations of motion display the fermionic pairing, 
i.e., bosonization, 

= i^Pn^, r = --^^pPi^^p. (155) 

Equation (|152p is the starting point for mean-field theory. The fermionic in- 
tegral is Gaufiian now. 



I- 



V(j)Vil}Vil}e-^^^ = I V4>e-^^, (156) 
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resulting in the purely bosonic action 

Sb-^ J d^x {m^(f>*(l) - lndet[i^ + h{PL<l) - Pr(I>*)]} ■ (157) 

Mean-field theory now neglects bosonic fluctuations and assumes that the 
bosonic ground state corresponds to that of the classical bosonic action. Of 
course, the Indet is still a complicated nonlinear and nonlocal expression; 
nevertheless, assuming that the ground state is homogeneous in space and 
time, (f) =const., the determinant can be computed by standard means [23j . 
For our purposes, it suffices to know that for 

A > ^ = Acr, (158) 

with A being the UV cutoff, the resulting effective potential Vb{4'*4') has a 
nonzero minimum, implying a nonzero vacuum expectation value (0) 7^ 0. In 
the fcrmionic language, this vacuum expectation value corresponds to a bi- 
fermionic condensate (</>) ~ (ipip) (a chiral condensate in the QCD context). 
The expectation value generates fermion mass terms ^ mt'ipj^ip witlF^ 

mi - h{(j)), (159) 

and the U(l)xU(l) symmetry is spontaneously broken to U(l) (fermion num- 
ber conservation); this implies the existence of one Goldstone boson corre- 
sponding to excitations of the phase of the (j) field. In the QCD context, this 
scenario corresponds to the spontaneous break-down of chiral symmetry with 
the pseudo-scalar mesons as Goldstone bosons. For A < Acr, the vacuum ex- 
pectation value is zero, ((/)) = and the system remains in the symmetric 
phase. 

These models of NJL type with broken symmetry show many similarities 
with QCD phenomenology, but there are important caveats from a microscopic 
viewpoint. First, there is no microscopic four-quark (or higher) self-interaction 
beyond criticality in QCD, A|yi 0. In other words, A[('0'(/')^ — ('/'TsV')^] and 
other four- fermion operators are RG irrelevant; the Euclidean microscopic 
action for vanishing current quark masses is given by 

^QCD = J d^x + i^w) , (160) 

where — dfj_ + igr'^A'^^ denotes the gauge covariant derivative in the funda- 
mental representation, which is generated by the with tr[T"T''] = {1/2)6"-''. 
Of course, four-quark operators in the effective action are generated by gluon 
exchange from fluctuations described by box diagrams; see Fig. [TT] The re- 
sulting (3 function for the four-quark coupling reads to lowest order 



The occurrence of 75 in the fermion mass term arises from our fermion conventions 
[100] : these are related to more standard conventions by a discrete chiral rotation. 
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Fig. 11. Box diagrams with fundamental QCD interactions generate effective four- 
fermion self-interactions interactions A. (Only one diagram per topology is shown; 
further diagrams with the regulator insertion (filled box) attached to other internal 
lines, of course, also exist.) The resulting flow-equation contribution to the running 
of A is given in Eq. (|161|l . 

dt\ = P\ = -cx ^ g", (161) 

where c\ is a coefficient which depends on the algebraic structure of the 
theory and the details of the IR regularization. For definiteness, let us con- 
sider QCD with an SU(3) gauge sector but only one massless quark flavor 
A^f = 1, the classical action of which has the same "chiral" symmetry proper- 
ties as the NJL model used above. For this system, the coefficient c\ obtains 
c\ = 5/(127r^) > for the regulator ([33| and a Fierz decomposition as chosen 
in [lOlj . Obviously, the four-quark self-interaction A is asymptotically free as 
it should be for a QCD-like theory. Of course, for increasing gauge coupling g 
towards the IR, a naive extrapolation of Eq. (|16ip predicts that the fermionic 
self- interaction can become critical A > Acr for some IR scale kcr- If this holds 
also in the full theory, the quark self-interaction X[{-ipip)'^ — (V'TsV')^] becomes 
strongly RG relevant at scales below fccr and we expect the system to end up 
in the symmetry-broken phase [102j . 

This scenario appears to match our expectations for QCD. But in order 
to arrive at a quantitative description we do not only face the problem of 
computational control in the nonperturbative domain; moreover, we have to 
deal with the conceptual problem of how to switch from the description in 
terms of quarks and gluons to that involving boson fields as well. In view 
of the Hubbard-Stratonovich transformation, we may be tempted to apply 
partial bosonization at some scale /cb < fccr- However, it turns out that this 
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leads to a strong spurious dependence on the precise choice of in generic 
truncations. This has to do with the following observation. Let us naively 
partially bosonize at < k^ . Diagrammatically, 




Sf =^ Sfb, (162) 

where A = in Sfb after bosonization. Now, let us perform another RG step 
and integrate out another momentum shell Ak. Owing to the box diagrams, 
new quark self- interactions are generated again in this RG step, 



k-Ak: dtX 



(163) 



In other words, the bosonizing field (f> which did a perfect job at k^ is no longer 
perfect at fc — Ak; it does no longer bosonize all quark self-interactions which 
it was introduced for. Incidentally, this problem does not only occur if gauge 
interactions are present, as Eq. (|163p seems to suggest. The same problem 
arises, e.g., in purely fermionic systems where fermion self- interactions are 
generated by the flow in many different channels also with nonlocal structure; 
since partial bosonization with a local bosonic interaction can never account 
for all four-fermion vertices, the remaining four-fermion structure will again 
generate the full structure in the RG flow. 

In a regime where many couplings run fast, neglecting the newly generated 
terms ()163|) can introduce large errors. But keeping these terms seems to make 
partial bosonization redundant. The solution of this dilemma is the subject 
of the following subsection. 

5.2 Scale-dependent field transformations 

In Eqs. (|162p . (|163p . we have observed that different bosonizing fields are 
needed to compensate the quark self-interactions at different scales. This 
points already to a solution of the problem [94j : we promote the bosoniz- 
ing field (f) to be scale dependent, — > ^j., the flow of which can be written 
as 

=Cfc[0,V,V5,...]. (164) 

Here, Ck is an a priori arbitrary functional of possibly all fields in the system. 
For the present problem, the idea is to choose Ck such that the resulting 
effective action does not possess fermionic self- interactions. For more 
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general cases, the functional can be chosen such that the effective action 
■ffe [0fe] becomes simple, since its simplicity is a strong criterion for the proper 
choice of the relevant degrees of freedom. 

We can formulate this idea in a differential fashion: we are looking for a 
functional Ck which yields a flow of /fc[(/)/c] taken at fixed (j)k, (suppressing 
further field dependencies on ^p^ip, . . 

dtFkiMU. = dtPkiM - I ^^^^ dtcf>k, (165) 



such that the flow of the fermion self-interaction vanishes for all fc, 9(A|0^ = 0. 
Therefore, if A = holds at one scale k, X stays zero at all scales, and (j)k 
becomes the "perfect" boson an all scales. The right-hand side of Eq. (|165p 
can be read as follows: the first term denotes the full RG fiow given in terms 
of a flow equation, whereas the second term with dt(f>k ~ Ck characterizes how 
(jik has to be modifled scale by scale in order to partially bosonizc fermion 
self-interactions on all scales. This fixes the functional form of Ck- 

Now, we need a flow equation for the effective action A [</>*:] with scale- 
dependent field variables. This can indeed be formulated in various ways. Here, 
we follow a general and flexible exact construction given in [8]. Consider the 
modified generating functional 

Zk[J] =e'^^[^l = / 2?^e-^['^l-^•/'^'=■^'''^''+•/'■^'^^ (166) 



where we have coupled a scale-dependent ipk to the source and the regulator. 
This combination guarantees that the resulting fiow equation has a one-loop 
structure p9l . The scale dependence of (p is given by 

dtipk ^ Ck[ip], (167) 

similar to Eq. p64p with the difference that Eq. (|167p is formulated under 
the functional integral, whereas Eq. (|164p holds for the fields conjugate to the 
source J, 

0, = {<f>k) ^ (168) 

Hence, Ck and Ck generally are not identical. 

The derivation of the flow of W^fe[J] is straightforward: 

dtWk[J] = -^^^ J '^f (^Jdtfk - ^ J 'PkdfRkfk - J V'kRkdtipk 



1 f s 

J{dt(pk) - -TidtRkGk - / jjRk{dtVk) 



(f>kRk{dtfk) - I (f>kdtRk(f>k- (169) 
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Here, we have defined the propagator similar to Eq. ([T 

Gkip) - (^) (P) = {^khp)^kip)) - M-p)Mp)- (170) 

We have also used the relation {(pkdtfk) = {jj + (t>k){dt'^k)- As usual, we 
define the effective action by means of a Legendre transformation, this time 
involving the scale-dependent field variables (cf. Eq. ((20)l ) , 

FkiA] = sup (^J Jcpk - Wk[J]^ - i y" <^fci?fc0fc. (171) 
The resulting flow of this effective action is 

dtFklM = \TYdtRkGk + j (Cfc^) Rk{dtVk) + j J^(9tc^k - {dt^k)), 

(172) 

where, in the course of the Legendre transformation, all J dependence turns 
into a cf)k dependence by virtue of J = Jgup — J[(f)k\') this also implies jj — 

For a given scale-dependent field transformation dti^k = C^Vp], we can 
successively work out {dt(pk), 4>k — (Vk) and dt(j)k = dt{^k)', in general, the 
latter is not identical to {dt(pk)- Following this strategy, we would only then 
be able to compute the flow of Fk [(j)k] ■ Also Ck would be a derived quantity, 
fixed implicitly by the choice of Ck ■ 

By contrast, we can supplement the flow equation ()172p with a bootstrap 
argument: since all we want to choose in the end is dt(l)k — Cfc [</>], the precise 
form of dtip ~ Ck need not be known; in fact, Ck does not occur directly in 
Eq. (|172p . but only in expectation values. Therefore, we simply assume that 
a suitable Ck exists for a desired Ck such that 

(dtipk) = dt(bk. (173) 

Of course, this is a highly implicit construction, and in view of the compli- 
cated structure of the mapping Ck — > Ck , the existence of a suitable Ck for an 
arbitrary Ck is generally not guaranteed or, at least, difficult to prove. Nev- 
ertheless, since the resulting flow equation will, in practice, be used together 
with a truncation, it is reasonable to assume that Eq. (|174p can at least be 
satisfied to the order of the truncation. As a consequence of Eq. (|174p . the 
flow equation simplifies, 

d,rk[A * d^FkiM - I ^ d,4>k (174) 

J 0(t)k 

^ Tvd.RkGk + [ ( Gk^) Rkdtcbk ~ [ ^-^^dt^k, 



2 



which is the desired flow equation for scale-dependent field variables. Apart 
from the standard first term ~ TrdfRkGk and the third term which carries 
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the explicit scale dependence of 0^, we encounter the second term which takes 
care of fluctuation contributions to the renormalization flow of the operator 
insertion dt^Pk in the functional integral. Actually, this second term will gen- 
erally be subdominant for not too large coupling: first, it is of higher order in 
the coupling, and second, Rk insertions lead to weaker numerical coefficients 
than dtRk insertions for standard regulators. We expect that this term does 
not induce strong modifications for couplings up to 0(1). 

5.3 Scale-dependent field transformations for QCD: 
Rebosonization 

Let us now turn back to our original problem of QCD-like systems, where 
quark self-interactions are generated by gluon exchange. In order to arrive at 
a mesonic description in the infrared, we now want to apply scale-dependent 
field transformations that translate the quark self-interactions into the bosonic 
sector an all scales - a process that may be termed partial rebosonization. 

In order to illustrate the formalism, let us study one-flavor QCD in a 
simple truncation. Apart from the standard kinetic terms for the quark and 
gluons, supplemented by wave function renormalization factors, we include a 
point-like four-quark self-interaction in the scalar-pseudo-scalar sectoJ"] 

r^.k = J d^x + iZ^^p^ + \ XUi^f - (^75^^)']) . (175) 

The initial condition of the four-quark operator in the UV is obviously given 
by A|i;_»yi — > 0. As a first step towards rebosonization, we include a complex 
mesonic scalar field in the truncation on equal footing, 

Fk = rp,k + J d'^x {z^d^<f>*d^<f> + y(0*0) + h[ijPi^<f>ib - v^Pr^V]) ,(176) 

with a scalar potential ^(0*0) = m'^(f)*(l) + 0{{(j)*(j)Y). The initial conditions 
for the scalar field at fc ^ yl need to be chosen such that the scalar has no 
observable effect on the QCD sector whatsoever. This is easily done by de- 
manding that the Yukawa interaction with the quarks vanishes h\k—,A —>- 0. 
We also choose a large scalar mass m^lk^A = 0{A'^), ensuring a fast de- 
coupling of the scalar; finally, we set Z^\k^A which makes the scalar 
non-dynamical at the UV scale. Solving the fiow with these initial conditions. 



Of course, in order to avoid any ambiguity with respect to possible Fierz rear- 
rangements of the four-fermion interactions in the point-like limit, all possible 
hnearly-independent four-fermion interactions, in principle, have to be included 
in the truncation. For simplicity, we confine ourselves here just to the scalar- 
pseudo-scalar channel, where chiral condensation is expected to occur. For the 
four-fermion interactions that will be generated by the flow, we use the Fierz 
decomposition as proposed in [1U1| . 
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the scalars rapidly decouple and only the standard QCD flow remains, reveal- 
ing the purely formal character of this first step towards rebosonization. 

As a second step, we now use the freedom to perform scale-dependent field 
transformations, as suggested in Eq. (|164p . We promote the field (f) to be scale 
dependent, and choose the functional Ck characterizing this scale dependence 
to be of the form 

dt(t)k = -ipPRil^dtak, dt(f>l = dtak, (177) 

with some function to be determined below. At this stage, let us study 
the consequences of the last term in Eq. ()174p on the resulting flow; note that 
we are now dealing with a complex field, such that this term goes over into, 
j ^dt(l)k ^ I ^dt(l)k + J ^dt(f>l. From the Yukawa interaction, we get 
from this term, together with Eq. (|177p . 

/i[V;Pl0V' - V^-Pr0>] ^ hdtakii^^p)^ - i^-/5^)% (178) 

which is a contribution to the flow of the quark self-interaction. Together with 
contributions from the first term of Eq. (I174p , and neglecting the second term 
of Eq. (I174p here and in the following as discussed above, we obtain the flow 
of the four-quark coupling at fixed 0^ , 

dt^Ut,^ -cx-^g'^ -hdtak, (179) 

where we used the result given in Eq. p6ip . Choosing the transformation 
function 

5ta, = -c.^^, (180) 

we obtain 

dtXW = 0, (181) 

which, together with the initial condition A|fc_,/i implies that A = holds 
for all scales k. The scale-dependent transformation Eq. (|177p thus has re- 
moved the point-like four-quark interaction in the scalar-pseudo-scalar sector 
completely by partial rebosonization. The information about this interaction 
is transformed into the scalar sector; for instance, the scalar mass term is also 
subject to the last term of Eq. ()174p : 

m^^*^ ^ ~m^dtak[4!Pi^H - V'^'r0*V']- (182) 

This yields a contribution to the flow of the Yukawa coupling. Together 
with the contribution from the first term of Eq. ()174p . i.e., the standard flow 
term, we obtain 

dth\4,^ = -^^Chg^h + w? dtak 
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(a) ^ (b) 

Fig. 12. Diagrams contributing to the RG flow of the scalar sector; (a) flow of the 
Yukawa coupling, see Eq. (|183|) : (b) flow of the scalar mass, see Eq. (|184|) . Only 
one diagram per topology is shown; further diagrams exhibit the regulator insertion 
(filled box) attached to other internal lines. 



where the coefficient is a result of the diagram shown in Fig. [T2[ a). For the 
linear regulator and in the Landau gauge, the result is Ch — I/tt^ for SU(3). 
The second term of Eq. p83p accounts for the fact that the fermions couple to 
a scale-dependent boson. The flow of the scalar mass term is not transformed 
by the choice of Eq. (|177p : only the standard flow-equation term contributes, 

dtm^ = Cmk^ h^, (184) 

where the coefficient c,„ yields for the regulator (|33p c„i — Nc/{8tt'^), resulting 
from the diagram in Fig. [T^b). The physical properties of the resulting bo- 
son field can best be illustrated with the convenient dimensionless composite 
coupling 



(185) 



and its /3 function 



dte = ~2e + c,„ + Chg'e + 2cxg^e^, (186) 



resulting from Eqs. (|16ip . (|183p . (|184p . The last term comes directly from re- 
bosonization. Since all Ci > 0, this /3 function looks like a parabola; see Fig.fT^ 
solid line. Without rebosonization, this (3 function would have corresponded 
to a straight fine, see Fig. [131 dashed line. 

The (3 function Eq. p86p exhibits two fixed points: e*_i is IR repulsive and 
e*^2 is IR attractive. In a small-gauge-coupling expansion, the positions of the 
fixed points are given by 



eV^^+0(5'), ~e.,2^^ + Oil/g'). (187) 
2 cxg^ 



Without rebosonization, only e*^i is present. 

If we start with initial conditions such that e\k^A < e quickly becomes 
negative, corresponding to the bosonic mass term dropping below zero, < 
0. This indicates that the potential develops a nonzero minimum, giving rise 
to chiral symmetry breaking and quark mass generation. However, we obtain 
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Fig. 13. Schematic plot of tlie /3 function (|186|) for tlie composite coupling e = -j^rj^ 
with arrows pointing along the flow towards the IR. The fixed point e,,i is IR 
repulsive; in its vicinity, the scalar field behaves as a fundamental scalar (dashed 
line). If the fiow is initiated with e|fc=A < ^ drops quickly below zero and the 
system runs into the regime with chiral symmetry breaking (xSB). For £|fc=A > 
the system rapidly approaches the bound-state IR fixed point e,,2, where the scalar 
exhibits bound-state behavior. QCD initial conditions correspond to e\k-tA oo. 



this initial condition ik^A < only if cither the scalar mass is small or the 
Yukawa coupling is large or both; but this is in conflict with our QCD initial 
conditions, specified below Eq. ()176p . With this initial condition, the system 
is not in the QCD universality class. Near e^^^i, the slope of the /3 function 
(|186p is —2, which is nothing but a typical quadratic renormalization of a 
bosonic mass term; the scalar behaves like an ordinary fundamental scalar 
here. In fact, there is a correspondence between e*^i and the critical coupling 
of NJL-like systems, 

Nc 

cf. Eqs. (|153p . (|158p . (|185p : the factor of A'c takes care of the additional color 
degree of freedom of the quarks which was not present in the NJL system of 
Subsect. 15.11 The initial condition e\k-tA < agrees with X\k^A > Acr, and 
the system is in the broken phase of the NJL model. 

For e\k^A > the system quickly approaches 2 either from above 
or below rather independently of the initial values of m,h\f^^A- QCD initial 
conditions with large initial m and small initial h correspond to e|fc^yi 00. 
But also for much smaller initial e, the system rapidly flows to e*^2, and the 
memory of the precise initial values gets lost. There, the system is solely 
determined by the gauge coupling which governs the fixed-point position. 
This is precisely how it should be in QCD. 
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Near e»_2, the boson is not really a fully developed degree of freedom. The 
flow does not at all remind us of the flow of a fundamental scalar, but points 
to the composite nature of the scalar. This justifies to call e*.2 the bound-state 
fixed point; for instance, in weakly coupled systems such as QED, the boson 
at the bound-state fixed point describes a positronium-like bound state. 

It is a particular strength of the RG approach with scale-dependent field 
transformations that one and the same field can describe bound state forma- 
tion on the one hand and condensate formation as well as meson excitations 
on the other hand; whether the field behaves as a bound state or as a funda- 
mental scalar is solely governed by the dynamics and the coupling strength of 
the system. 

So far, our analysis of the system was essentially based on weak-coupling 
arguments, revealing that QCD at initial stages of the flow approaches the 
bound-state fixed point. But, we still have to answer a crucial question: how 
does QCD leave the bound-state fixed point and ultimately approach the 
chiral-symmetry broken regime? The answer is again given by the gauge 
coupling which controls the whole flow. For increasing gauge coupling, the 
parabola characterizing the (3 function (|186p for i is lifted, as depicted in 
Fig. [HI At a critical coupling value, — g^r, the fixed points e*,i and e*,2 an- 
nihilate each other and the system runs towards the chiral-symmetry broken 
regime. This transition is unambiguously triggered by gluonic interactions. 
For instance in the present case of one-flavor QCD with gauge group SU(3), 

2 

the critical couphng is given by aa = ~ 0.74 for the regulator Eq. ([33]) . 
Since this coupling value is not a universal quantity, one should not overem- 
phasize its meaning; however, it is interesting to observe that this coupling 
strength is in the nonperturbative domain, as expected, but not very deeply. 
In particular, since loop expansions go along with the expansion parameter 
q/tt, the critical coupling and thus the approach to chiral symmetry breaking 
appears still to be in reach of weak-coupling methods (not to be confused with 
perturbation theory). 

We conclude that the continuous scale-dependent translation allows for a 
controllable transition between microscopic to macroscopic degrees of freedom 
and between different dynamical regimes of a system. From a quantitative 
viewpoint, no spurious dependence on a bosonization scale, i.e., a scale at 
which degrees of freedom are discretely changed, is introduced, because field 
transformations are continuously performed on all scales. This helps main- 
taining the predictive power of truncated RG flows. As a result, macroscopic 
parameters can quantitatively be related to microscopic input. 

5.4 Fiirther reading: aspects of field transformations* 

The scale-dependent field transformations introduced above were illustrated 
with the aid of rt-point interactions in the point-like limit, e.g., a four-fermion 
interaction in the zero-momentum limit. Of course, the formalism can also be 
used if momentum dependencies of the vertices are taken into account. This 
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Fig. 14. Schematic plot of the /3 function (|186|) for the composite coupling e — -j^fj^ 
with arrows pointing along the flow towards the IR. At weak gauge coupling, QCD- 
like systems first flow to the bound-state fixed point e,^2 where they remain over 
a wide range of scales. For increasing gauge coupling g^, the /3 function is lifted 
(dashed lines). At the critical coupling = g^^, the fixed points are destabilized 
and the system rapidly runs into the chiral symmetry broken regime (xSB). 

is particularly important in the context of rebosonization, since composite 
bosonic fluctuations and bound states manifest themselves by a characteristic 
momentum dependence in the fermionic n-point correlators; for instance, a 
bosonic bound state corresponds to a pole in the s channel of the fermionic 
Minkowskian 4-point vertex. 

A generic momentum dependence of an n-point vertex is nonlocal in co- 
ordinate space. By means of a field transformation, certain nonlocalities can 
be mapped onto a local description. The Hubbard- Stratonovich transforma- 
tion is an example for such a transformation which maps a specific nonlocal 
four-fermion vertex onto a local bosonic theory with a local Yukawa coupling 
to the fermion. 

With the following scale-dependent field transformation, this nonlocal-to- 
local mapping can be performed on all scales for a four-fermion vertex with s 
channel momentum dependencies, A — A(s — q^), 

dtMq) = i^PR^Mdtakiq), dtrkiq) = ~i^pPL^K-q)dtak{q).{189) 

This results in a flow of the 4-point interaction for the transformed fields 
which reads 

dtXU.^dtX-hdtakiq). (190) 

Obviously, the q dependence of ak can be chosen such the fluctuation-induced 
s-channel q dependence of A is eaten up for all scales and all values of q. 
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Via the momentum-dependent analog of Eq. (|183p , this induces a momentum 
dependence of the Yukawa interaction hb h{q); this can be taken care of 
by a further generalization of the field transformation, 

dtclykiq) = {i'PR^)iq) dtUkiq) - (3k{q)<t>k{q), 

dt(bUq) = -(V;PL^)(-g) dtakiq) - PkiqWkiq)- (191) 

The resulting flows for the Yukawa coupling and the inverse scalar propagator 
for the transformed fields are then given by 

dth{q)\^, = dth{q) + ^^f±^dtXiq') + hdMq), 
{dtZ^{q)q^ + dtm')\^^ = dtm" + 2dtl3k{q) [Z^q^ + m^). (192) 

For a given momentum dependence of A, the first equation can be used to 
determine the transformati<;n function Mq)B The second equation then fixes 
the running of the mass and of the wave function renormalization of the 
transformed field. In this manner, specific nonlocal momentum structures can 
be transformed from the fermionic interactions into a local boson sector. 

This strategy has been used for the abelian gauged N JL model [94 , QED 
[91] and QCD-like systems in [lOlj . At weak coupling, one typically finds that 
the bound-state fixed point discussed above of the composite coupling e has 
similar counterparts in many other couplings, such as the scalar dimensionless 
mass m? /k"^ and the Yukawa coupling h. This is a manifestation of their RG 
irrelevance at weak coupling, with the dynamics of the scalar sector being fully 
controlled by the fermions and their gauge interactions. If the gauge coupling 
becomes large, the bound-state fixed points in all these couplings is desta- 
bilized and interactions involving composite bosons can become relevant or 
even dominant close to a transition into a broken-symmetry regime. Since this 
switching from irrelevant to relevant is a smooth process under the flow being 
controlled by the dynamics itself, the predictive power of the computation is 
maintained. 

In QCD, the resulting effective action at IR scales naturally exhibits a 
sector which is similar to a chiral quark-meson model [103|, [81] I104j but with 
all parameters fixed by the outcome of the RG flow with field transformations. 
Also the gluonic sector can still be dynamically active and contribute further 
to the running of the chiral sector. 

A further application of the scale-dependent field transformation can be 
found in |105j where it is shown that the Fierz ambiguities mentioned in the 
preceding section can be overcome by treating all possible interaction channels 
and the corresponding scale-dependent bosonic composites on equal footing. 



The momentum-independent part can, for instance, be fixed such that dtZ^{q — 
k) = 0, ensuring that the approximation of a momentum-independent Z^ is self- 
consistent. 
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